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ABSTRACT 

A  laboratory  study  has  been  made  o-f  the  characteristics 
of  an  interf erometric  type  optical  fiber  pressure  gradient 
hydrophone.  The  optical  system  is  configured  as  an 
all -fiber  Mach-Zehnder  interferometer  excited  by  a  single 
mode  gallium  arsenide  diode  laser.  A  pair  of  identical 
fiber  coils,  one  in  each  arm  of  the  interferometer,  forms 
the  sensing  portion  of  the  gradient  hydrophone.  Each  coil 
consists  of  10  meters  of  polyethylene  jacketed,  single  mode 
fiber,  wound  in  a  doughnut  shaped  element  of  mean  diameter  4 
cm  and  thickness  3  mm.  An  acoustic  calibrator,  similar  in 
design  to  the  USRO  type  G40  shipboard  calibrator  used  by  the 
U.S.  Navy,  was  constructed  for  use  in  obtaining  single  and 
dual  (dipole)  coil  sensitivities.  The  acoustic  sensitivity 
of  each  coil  was  determined  separately  and  then  the  pair  was 
aligned  coax i ally  and  separated  by  10  cm  to  operate  as  a 
pressure  gradient  device.  Details  of  the  construction  of 
the  system,  calibration  procedure,  lock-in  demodulator,  and 
hydrophone  sensitivity  data  are  presented.  The  system  has 
demonstrated  a  14  dB  sensitivity  advantage  over  the  USRD 
type  G63  standard  (piezoelectric)  pressure  gradient 
hydrophone. 
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I ■   INIBQDUCT I ON 
A.   BACKGROUND 

The  use  of  optical  -fiber*  -for  communication  and  sensing 
is  a  rapidly  expanding  area  of  recent  technology.  Optical 
-fiber  technology  is  a  little  over  a  decade  old,  starting 
with  the  first  demonstrations  of  low-loss  (less  than  1000 
dB/km)  optical  fiber  waveguides.  Technological  advances  in 
material  processing,  fabrication  of  optical  fibers, 
interconnection  devices,  cables,  sources,  and  detectors  now 
allow  repeater less  light  transmission  over  great  distances 
(tens  of  kilometers)  with  an  extremely  large  bandwidth  CRef . 
13. 

Capitalizing  on  the  availability  of  optical  fiber 
components,  there  has  been  significant  interest  and  progress 
in  the  past  half  decade  in  exploitation  of  optical  fibers  as 
transduction  elements  in  a  new  class  of  sensors.  These 
sensors  ars  capable  of  detecting  acoustic  fields,  linear  and 
rotational  acceleration,  temperature,  strain,  electric  and 
magnetic  fields,  and  other  physical  parameters.  This 
versatility  can  also  lead  to  difficulties  and  much  effort 
has  been  expended  to  enhance  sensitivity  to  the  parameter  of 
interest  while  desensitizing  the  fiber  to  all  other 
fluctuations   CRef.  21.       Fiber   sensors  offer  a   number   of 
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potential  advantages:  increased  sensitivity,  geometric 
versatility  as  -Fiber  sensors  can  be  con-figured  in  arbitrary 
shapes,  lightweight  components,  and  low  cost.  Progress  in 
demonstrating  these  advantages  has  been  substantial  in  the 
past  -Few  years  with  over  60  di-f-ferent  types  of  sensors  using 
optical  -Fiber  as  the  transduction  element  being  developed. 
These  sensors  range  -From  simple  on/o-F-F  -Fluid  level 
indicators,  to  the  more  sophisticated  inter-Ferometric 
con-Figurations.  This  large  number  o-F  individual  devices  is 
usually  categorized  into  amplitude  or  phase  sensors.  In 
e-F-Fect,  the  sensor  modulates  some  parameter  o-F  the  lightwave 
in  an  optical  fiber,  such  as  the  intensity  or  the  phase.  In 
the  former  case  the  physical  perturbation  interacts  with  the 
fiber  to  directly  modulate  the  intensity  of  the  light  in  the 
fiber.  In  the  latter  case  the  physical  perturbation 
modulates  the  optical  phase  of  coherent  light  propagating 
through  the  fiber.  Optical  interf erometry  must  be  utilized 
to  convert  the  optical  phase  modulation  to  optical  intensity 
modulation  in  order  to  permit  photodetection.  This 
technique  allows  for  extremely  high  sensitivity.  Phase 
modulation,  conversion  of  phase  modulation  to  intensity  by 
interf erometry,  and  subsequent  detection  by  photodiodes  will 
be  discussed  in  more  detail  in  Chapter  II. 
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B.   PURPOSE  OF  THIS  STUDY 

Fiber  optic  acoustic  sensors  ware  -first  demonstrated  to 
be  feasible  in  1977  CRef.  3  and  43.  Since  then  a  significant 
degree  o-f  progress  has  been  achieved,  both  in  the 
understanding  of  the  relevant  acousto-optic  transduction 
mechanisms  as  well  as  in  component  development  and  sensor 
packaging  CRef.  5  and  63.  Most  of  the  fiber  optic  hydrophone 
work  performed  to  date  has  been  with  monopolar 
omnidirectional  devices.  Figure  1.1  illustrates  the  basic 
geometry  of  a  typical  fiber  optic  hydrophone  which  utilizes 
acoustically  induced  phase  modulation  for  transduction.  A 
laser  beam  is  split  into  two  arms  of  an  interferometer  where 
a  sensing  coil  in  one  arm  is  exposed  to  the  acoustic 
pressure  field  while  a  reference  coil,  forming  the  second 
arm,  is  isolated  from  the  field.  The  two  beams  are 
recombined  and  allowed  to  interfere  on  the  surface  of  a 
photodetector . 

Since  their  invention  in  the  late  1970s  two  unique 
characteristics  of  interf erometric  type  fiber  optic 
hydrophones,  i.e.,  their  basic  two-arm  design  and  their 
proven  high  sensitivity,  have  suggested  their  use  as 
gradient-type  devices.  A  fiber  optic  gradient  hydrophone 
has  a  sensing  coil  in  each  arm  of  the  interferometer  and 
both  are  exposed  to  the  pressure  field.  Figure  1.2 
illustrates  the  basic  geometry   of   an   acoustic  fiber  optic 
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gradient  hydrophone.  It  was  the  aim  of  this  project  to 
construct  and  partially  evaluate  in  a  laboratory  environment 
the  characteristics  of  a  fiber  optic  gradient  hydrophone. 

C.   FORMAT  OF  THE  REPORT 

The  theoretical  bases  of  conventional  gradient 
hydrophones,  the  calibration  of  gradient  hydrophones,  and 
the  behavior  of  interf erometric  sensors  used  in  preparation 
of  this  report  is  presented  in  Chapter  II.  Details  of  the 
construction  of  the  optical  fiber  sensor  »rm  presented  in 
Chapter  III.  Included  in  Chapter  III  is  a  description  of  the 
characteristics  of  the  individual  parts  of  the 
interferometer.  The  parts  are  discussed  in  the  order  that 
they  were  built  into  the  system.  Specific  instrumentation 
used  to  obtain  data  is  also  discussed  in  Chapter  III.  The 
experimental  procedures  used  to  establish  the 
characteristics  of  the  system  and  conduct  sensitivity 
measurements  of  the  individual  sensor  coils  and  the  gradient 
(dual  coil)  hydrophone  are  discussed  in  Chapter  IV.  Analysis 
of  the  data  and  evaluation  of  the  results  »rm  also  presented 
in  Chapter  IV.  Chapter  V  contains  concluding  remarks  and 
recommendations  for  further  work.  Appendix  A  lists  a  BASIC 
computer  language  program  used  in  data  acquisition  and 
Appendix  B  presents  an  example  of  the  data  obtained  using 
this  program. 
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I I -   THEORY 
A.   CONVENTIONAL  GRADIENT  HYDROPHONE 

In  many  acoustic  applications  it  is  necessary  to 
determine  the  direction  of  signal  arrival  as  well  as  to 
detect  acoustic  pressure  levels.  Directional  -Features  »ra 
usually  achieved  with  the  use  o-f  multiple  hydrophones 
spatially  distributed  in  some  specified  orientation,  e.g.,  a 
line  array.  The  simplest  o-f  these  designs  is  the  dipole 
pressure  gradient  hydrophone. 

Some  pressure  gradient  hydrophones  Ar9  designed  so  that 
the  sensing  part  o-f  the  hydrophone  oscillates  in  the  medium 
under  the  influence  of  the  sound  pressure  gradient.  That 
is,  a  hydrophone  is  constructed  for  which  the  driving  force 
and  the  hydrophone's  output  is  proportional  to  the 
difference  between  the  pressures  at  two  different  locations 
in  an  acoustic  sound  field.  Therefore,  the  electrical 
output  is  proportional  to  the  pressure  gradient  in  the  sound 
field.  The  main  characteristic  and  advantage  of  a  pressure 
gradient  hydrophone  is  the  directivity  pattern  obtained  from 
such  a  small  .size  hydrophone.  As  will  be  shown  below, 
pressure  gradient  hydrophones  have  dipole,  or  figure-eight, 
directivity  patterns,  as  indicated  in  Figure  2.1;  hence  they 
are  bidirectional.   Assuming   the   size  of  the  hydrophone  is 
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Figure  2.1     Directivity  Pattern  of  Pressure 
Gradient  Hydrophone 
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small  compared  to  the  acoustic  wavelength  \  ,  the  dipole 
response  when  oriented  at  any  angle  Q  relative  to  an 
incoming  plane  pressure  wave  is  proportional  to  cos Q . 

Both  for  simplicity  and  since  they  resemble  the  fiber 
optic  system  under  study,  pressure  gradient  hydrophones 
formed  by  combining  two  small,  closely  spaced,  individual, 
pressure  hydrophones  will  be  discussed.  The  performance  of 
this  hydrophone  pair  (or  doublet)  is  easily  calculated,  and 
the  behavior  thus  derived  applies  also  to  pressure  gradient 
hydrophones  of  more  complicated  construction  as  long  as  the 
assumptions  of  compactness  made  in  the  derivation  are  not 
violated. 

As  an  introduction  to  the  theory  of  this  type  of 
hydrophone,  consider  two  small  pressure  hydrophones  placed 
in  a  standing  acoustic  wave  field,  P(x,t),  separated  by  a 
distance  d  <<  A  as  shown  in  Figure  2.2.  The  hydrophones  are 
placed  so  that  a  null  in  the  standing  wave  pattern  occurs 
halfway  between  the  two  hydrophones.  The  electrical  outputs 
of  the  two  Arc  connected  180  degrees  out -of -phase  with  each 
other  to  obtain  a  dipole.  The  dimensions  of  the  two 
hydrophones  are  typically  much  less  than  a  wavelength, 
therefore  the  acoustic  field  is  assumed  to  be  negligibly 
disturbed. 
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Figure  2.2     Geometry  Used  in  Deriving  Sensing 
Characteristics  of  Acoustic  Dipole 
(Pressure  Gradient) 
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Consider  the   plane  sinusoidal  wave  shown  in  Figure  2.3. 

The  instantaneous  acoustic  pressure,  P(x,t)  is  given  by: 

jCJt 
P(x,t)  -  P  sinCkx3e 

o 

(2.1) 

where  P  is  the  peak  acoustic  pressure,  k  is  the  propagation 
wave  number  k  =  2  7T/A.*  and  x  is  the  distance  o-f  one  of  the 
hydrophones  from  absolute  minimum  pressure  (i.e.  the 
standing  wave  pressure  node).  Assuming  the  argument  o-f  the 
sine  -function  is  small,  the  -first  order  (linear) 
approximation  (i.e.,  sin  y  =  y)  o-f  the  pressure  amplitude  at 
x  is  given  by: 

jCJt 
P(x,t)  -  P  kx  e 

o 

(2.2) 

If  one  sensor  is  at  x  =  -d/2  and  the  other  at  x  *  +d/2  as 
shown  in  Figure  2.3,  then  the  amplitude  o-f  the  pressure 
difference  to  be  measured  may  be  expressed  as: 


Ap  -  P     -  P 

+d/2    -d/2 

(2.3) 
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Figure  2.3    Pressure  Distribution  as  Function  of 
Distance  from  Null  Pressure  Point 
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This  may  be  written  as: 


Ap  =  P  kCd/2  -  <-d/2)D 
o 


(2.4) 


Or: 


Ap  ■  P  kd 

o 


(2.5) 


For  the  more  general  one  dimensional  case  of   any   point 

along  a  standing  wave,  consider  a  standing  wave  in  a  tube  o-f 

water  as  shown  in  Figure   2.4.   The   instantaneous   acoustic 

pressure  P(x,t)  is  given  by: 

JCJt 
P(x,t)  =  P  sin(kx)e 

o 

(2.6) 

where   x   has  been   chosen   to   be   zero   at   the   air-water 
boundary. 

I-f  we  expand  equation  (2.6)  into  a  Taylor's  series  about 
x  ■  h,  where  h  is  the  distance  from  the  air-water  boundary, 
the  following  equation  results: 

(x,t)  =  P(h,t)  J 3  P(h,t)\(x  -  h)  Y(3  P(h,t)\(x  -  h)  +.  .  . 

£«      J  W    J~T~ 

(2.7) 
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Figure  2.4    Pr assure  and  Velocity  Distribution  in 
Standing  Wave  Tubs. 
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where  the  first  term  and  all  their  partial   derivatives   are 

evaluated  at  x  =  h.   Time  t  is  held  constant. 

The  -First  partial  derivative  (pressure  gradient)  iss 

jOJt 
ccas(kx)i 

o 


UP/dx    =  P  kcos(kx)e 


(2.8) 


The   second   partial    derivative   isi 

^2  2  2  jCJt 

0   P/Qx      -   -P   k   sin(kx)e 


(2.9) 


The   pressure   difference   between   two   points   ax i ally 
displaced  from  one  another  by  distance  Ax  where  /\  x  =  x  -  h 


is: 


Ap  =*  P(x,t)  -  P(h,t) 


(2.10) 


or 


Ap  -  (5p/6x)Ax  +  (d  p/(3 x  )(Ax)  /2  +  ... 


(2.11) 


If  l\x  is  small  compared  to  the  wavelength  of  the  standing 
wave  all  terms  above  second  order  may  be  neglected  and 
therefore  the  series  may  be  truncated  after  the  second  order 
term. 
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If  h  and  A  x  are  selected  so  that  a  pressure  node  is 
located  at  Ax/2,  then  the  difference  in  pressure  /\p 
reduces  to  equation  (2.5)  with  Ax  ■  d,  in  the  linear 
approximation.  The  first  correction,  6  ,  to  the  linear 
approximation,  in  the  vicinity  of  a  pressure  node,  becomes i 

6  -  -C<kAx>/23tan<kh> 

(2.12) 

A  maximum  in  the  pressure  difference  will  occur  when 
each  sensor  is  equally  distant  from  the  pressure  node  and 
conversely  a  minimum  in  the  pressure  difference  will  occur 
when  each  sensor  is  equally  distant  from  a  pressure  maximum 
( ant i node) . 

A  pressure  gradient  hydrophone  is  sometimes  called  a 
velocity  hydrophone,  the  implication  being  that  it  responds 
to  the  particle  velocity  in  the  acoustic  medium.  This  is 
somewhat  misleading  since  the  velocity  of  the  hydrophone 
surface  is  normally  much  smaller  than  the  particle  velocity 
of  the  incident  acoustic  wave.  But  even  though  the 
hydrophone  is  not  a  direct  detector  of  particle  velocity, 
its  response,  being  proportional  to  pressure  gradient,  is 
also  proportional  to  particle  velocity  of  the  incident 
acoustic  wave  as  shown  in  the  one  dimensional  linearized 
Eu 1 er  equat  i  on : 

-(3p/dx  -  JCjpu  (2.13) 
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where  (a)  is  the  angular  -frequency  of  the  acoustic  wave,  Q 
the  mass  density  of  the  medium,  and  u  is  the  Eulerian 
particle  velocity  of  the  free-field  acoustic  signal. 

B.   CALIBRATION  OF  GRADIENT  HYDROPHONES 

Pressure  gradient  hydrophones  may  be,  and  usually  ^rm, 
calibrated  in  terms  of  pressure  rather  than  velocity 
(pressure  gradient  is  never  used  as  a  reference).  Since  the 
relationship  between  pressure  and  velocity  depends  on  the 
wave  and  boundary  conditions,  these  conditions  must  be 
known,  be  standard,  or  fit  the  definition  of  sensitivity. 
Therefore  the  sensitivity  of  pressure  gradient  hydrophones 
usually  is  given  in  terms  of  volts/micropascals  (V/LtPa) 
specified  at  a  particular  frequency. 

Plane  progressive  waves  are  specified  in  the  definition 
of  free-field  voltage  sensitivity,  and  the  relationship 
between  the  pressure  p  and  velocity  u  in  such  waves  is  given 
by  the  wave  or  characteristic  impedance  p/u  =  He.  Thus,  the 
sensitivity  of  a  pressure  gradient  hydrophone  expressed  in 
terms  of  pressure  differs  from  the  sensitivity  expressed  in 
terms  of  velocity  by  the  constant  factor  Dc,  where  D  is  the 
density  of  the  medium  and  c  is  the  speed  of  sound. 

The  standard  procedure  for  pressure  gradient  hydrophone 
calibration  is  by  direct  comparison  with  a  standard  pressure 
hydrophone  only  if   essentially   plane  waves  impinge  on  both 
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hydrophones.    Since    waves   in   a   free   -field   never   are 
perfectly  plane  or  if  a  spherically  diverging   wave   from   a 

point     source     is     used,     a     correction     factor, 

\2       2 
101ogCl+  A  /(27Tr)  1,   where   r   is  the  distance  from  the 

source,   must   be  subtracted  from  the   measured   sensitivity 

CRef.  73.  This  correction  factor  decreases  as  r  increases. 

Because  of   the  difficulties   in   obtaining   free-field 

conditions  particularly  at  low  frequencies,  a   standing-wave 

tube  like  that  shown  in  Figure   2.4  was   employed   for   the 

laboratory  investigations  described  in  this  report.   Since  a 

water-air   boundary   is   a   near   zero   impedance  boundary, 

essentially  complete  reflection   will   take  place  at  the  top 

of  the  tube.   A  standing-wave  system  then  can  be  established 

and  the  following  relations  hold  in   the   ideal   case   (i.e. 

SWR=  OO  )  s 


p  =  p  sin  kh, 

o 

(2.14) 


u  -  (p  /Qc)cos  kh, 

o 

(2.15) 


p/u  =pc 


tan  kh. 


(2.16) 
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I-f  the  pressure  is  measured  at  a  point  a  distance  h  -from 
the  water-air  surface,  equations  (2.14)  and  (2.15)  can  be 
used  to  compute  the  pressure  and  velocity  at  any  other 
point.  As  stated  earlier  the  hydrophones  are  assumed  to 
have  a  negligible  effect  on  the  standing  wave  pattern.  If 
the  tube  is  operated  below  cutoff  for  non-axial  modes,  then 
the  pressure  variations  normal  to  the  tube  axis  can  be 
neglected.  A  detailed  discussion  of  the  calibrator  used  is 
given  in  Chapter  III  section  D7. 

A  more  detailed  explanation  of  sensitivity  and 
calibration  procedures  can  be  found  in  Bobber  CRef  73.  In 
this  experiment  sensitivities  of  single  coil  fiber  optic 
acoustic  hydrophones  are  defined  in  terms  of  fJi  radians///  Pa. 
Sensitivities  of  gradient  (dual  coil)  hydrophones  are 
defined  in  terms  of  optical  phase  shifts  induced  by  acoustic 
pressures  and  are  expressed  in  fl  radians/  fj. Pa/cm.  The 
procedures  to  obtain  fiber  optic  hydrophone  sensitivities 
(single  and  dual  coil  designs)  are  discussed  in  Chapter  IV. 
A  thorough  presentation  on  fiber  optic  acoustic  sensitivity 
is  given  by  Davis,  et  al . ,  CRef.  63. 

C.   FIBER  OPTIC  ACOUSTIC  SENSOR 

1.   Acoustically  Induced  Phase  Modulation 

Laser  light  conducted  by  optical  fibers  submerged  in 
a  liquid  medium  may  be  modulated  (intensity  or  phase)  by 
acoustic  pressures.   This  phenomenon  is  exploited  to  make  an 
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acoustooptic  sensor  of  underwater  sound  by  demodulating  the 
light  using  a  specific  detection  scheme.  Carome  and 
Satyshur  CRef.  83  present  in  detail  general  fiber  mode 
theory  and  the  theory  of  phase  and  intensity  modulation 
induced  in  optical  fibers.  Phase  modulation  occurs  when  the 
instantaneous  phase  angle  of  the  lightwave  in  the  optical 
fiber  is  varied  proportionally  in  accordance  with  the 
instantaneous  value  of  a  modulating  signal,  e.g.  a 
si nusoi dally  varying  acoustic  field. 

Two  major  detection  schemes  for  inter f erometric  type 
sensors  discussed  in  the  references  are  homodyning  and 
heterodyning  CRef.  5  and  63.  Homodyne  detection  in  optical 
interf erometry  implies  that  the  reference  beam  has  the  same 
optical  frequency  as  the  sensing  beam.  On  the  other  hand, 
heterodyne  detection  is  based  on  the  use  or  mixing  of  two  or 
more  frequencies,  i.e.,  the  frequency  of  the  reference  and 
sensing  beams  differ  by  l\(x)  .  Interf  erometry  is  used  in 
either  case  to  recover  the  acoustic  signal  frequency. 
Homodyne  detection  was  used  throughout  this  experiment  so 
only  it  will  be  considered  in  detail. 

An  external  pressure  field  applied  to  a  fiber  will 
cause  various  changes  in  the  fibers  physical 
characteristics,  such  as  changes  of  the  core  radius  and 
length,  and  changes  in  the  values  of  the  optical  indices  of 
refraction  of  the  core  and  cladding  CRef.  3  and  43.  For  the 
case  of  pressure  induced  changes  of  index  and  of  length,  the 
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phase  shift  A,  produced  by  a  change  of  pressure  A  P  might  be 
given  asi 

A  -  k  PJ^C<(3n/(3p>  +  n/J^<  (3jP/(3p>:1 

o 

(2.17) 

where  n  is  the  optical  index  of   refraction  of  the  core,  k 

o 

is  the  propagation  constant  of  light  in  the  fiber,  P  is  the 
acoustic  pressure,  and  y  is  the  length  of  fiber  subjected  to 
the  pressure.  The  partial  derivative  (Qn/CJP)  is  evaluated 
at  constant  entropy.  Whatever  the  form  of  the  bracketed 
term  of  equation  (2.17),  the  phase  shift  amplitude  A  enters 
the  analysis  in  identical  ways  so  that  it  is  possible  to 
proceed  without  being  concerned  with  its  specific  form. 

Using  a  single  frequency  coherent  laser  source,  the 
time  variation  of  the  electric  field  vector  of  the 
lightwave,  propagating  in  a  fiber  that  has  been  subjected  to 
an  acoustic  pressure  field,  may  be  expressed  asi 

E(t)  =  E  expCjCCJ  t  -i-  AsinC  CJ  t)3 
o         o  s 

(2.18) 

where  CJ   i*  the  angular  frequency   of   the   coherent   laser 

o 
source,  (jj    ,  the  angular  frequency  of  the  sound  field,  and  A 

s 
is  the  phase  shift  amplitude  as  described  above. 

Phase  modulation   cannot   be   directly   detected   by 

photodetectors  due  to  the  fact  that  they  cannot   follow   the 


31 


instantaneous   values   o-f   such   high   -frequency   variations 

14 
(approximately  10    Hz).  Thus,  in  order  to  accomplish  phase 

detection,  an   inter-f erometric   technique  must   be  used  to 

convert   the  phase  modulation   to   intensity    (amplitude) 

modulation   prior  to  detection.   A  monochromatic  laser   beam 

is  split,   with  part  o-f  the  beam  being  transmitted   into  a 

re-ference  -fiber.   The  other  beam   is   sent   into   a  sensing 

-fiber  exposed  to  an   acoustic   -field.    Pressures  exerted  on 

the  sensor   -fiber  as  a  result  o-f  the  sound   energy  density 

produce  optical  path  length  changes  in  the  fiber   which,   in 

turn,   cause  phase  modulations  in   the   lightwave.    A-fter 

passing   through   the   -fibers,    the   two   beams   are   then 

recombined.   Thus  the  total  electric  -field  is  given  byi 


E   -  E  (t)  +  E  (t) 
T     1        2 

(2.19) 


where 


E    (t)    -  E  expCj(CJ  t   +  Asin(CJ  t)3 
1  1  o  s 

(2.20) 


and 


E    (t>    *_EkexpCj((J  t   +0  )3 
2  2  o 


(2.21) 
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where  (7)  is  an  arbitrary  phase  difference   which   will   vary 

with   temperature,   the  optical  path  difference,   and   other 

factors.   E  (t)  is  the  electric  field  vector  of  the   sensing 

_>  1 
arm  and  E  <t>  is  the  electric  field  vector   of  the  reference 

2 
arm.    The   intensity    (I)    of   the  recombined   beams   is 

proportional   to  the   magnitude   of   the   square  of   E  . 

T 
Therefore,   after   mathematical   manipulation   as   shown   by 

Carome  and  Satyshur  CRef.  83  and  neglecting  terms  that   vary 

at  angular  frequency  GJ   and  2  (jj   *  since  as   stated   earlier 

o        o 
they  are   undetectable   by   a   photodetector ,   the  following 

equation  resultsi 

I  CXEj    +  E,  +  E  •  E  cos(Y)J  <A> 
2     2         2       0 

+   2E,«  E1sin0J(A)sinCL)t    +   2E/-  Ezcos0J    <A)cos2(Jt 

+   2e"«  *E   sinC^J,  (A)sin36Jt    +    .    .     .    . 
'       l  ~  3  S 

(2.22) 

where  E  •  E   makes  explicit  the  fact   that   the   polarization 

1   2 
directions  may  not  be  the   same.    The  first  three  terms  are 

considered  the   DC   terms   and  all  other  terms  of  the  series 

»ra   the  AC  terms.   Thus   the   resulting   intensity   function 

consists  of   a   series   of    harmonics   of   the   acoustic 

frequencies  where  the  amplitude   of  each  successive  harmonic 

is  a   function   of   the  acoustic  pressure,  and  varies  as  the 

Bessel  function  of  corresponding  order. 


33 


If  both  of  the  beams   are   phase   modulated   by  the 

acoustic  pressure  as  in  the  case  of   the   pressure  gradient 

hydrophone,  with  phase  shifts   of   A    and  A  ,  respectively, 

1        2 
the   final   expression   for   the  intensity  is   identical   to 

equation   (2.22),    except    that    for    A,    the   quantity 

A A  a  A  -  A   must  be  substituted. 
1    2 
The  recombined  lightwave  with  its  varying   intensity 

is   then   projected   onto   the  surface  of   a   photodetector. 

There  results  in  the  photodetector  current  a  signal  which  is 

related   to  the   acoustically   induced   phase   shift.    The 

photodetector  current  has  components  of  the  following  form: 

00 


i(t)    -   i    cos<j6*    iJ    (k   x)    +   2  Vj       (k   x)cosC2n  ((J  t>  3> 
0  0      o  n  —  i    2n      o  s 

oo  "-/ 

-i    sin0*    (2  Vj  (k   x)sinC (2n+l) (GJ t) 3> , 

0  n-r,  2n+l      o  s 


n=o 

(2.23) 

where  J   is  the   Bessel   function   of  order  n  CRef.  83,  <T)  is 

n 
the    static    phase   difference,   which   will    vary    with 

temperature,  path   length   difference,   and   other   factors, 

k  =2 7T/  \    ,  where   k   is  the  optical  propagation  constant  in 

o        o  o 

the  fiber,  \       is  the  optical  wavelength  in  the  fiber,  and  x 

o 
is  the  amplitude  of  the  optical  path-length   change   at  U)    - 

s 
The   J    term   in   equation   (2.23)   is   the   DC    component 

0 
associated  with   the  acoustically   induced  signal   with   the 
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remaining  terms  being  the  AC  component.   Only  the   AC   terms 

of    equation   (2.23)   will   be   considered   in    succeeding 

sections. 

2-   Feedback  Stabilization 

The  purpose  of  homodyne  detection  in  the  -fiber  optic 

interferometer   discussed  above  is  to  trans-form  the   optical 

output   of   the   interferometer  into  an   electrical   signal 

proportional  to  the  amplitude   of   the   relative  phase  shift 

produced   by   a   pressure   oscillation   at   frequency  U)    . 

s 
Homodyne  signal  processing   suffers   from  phase  drift  due  to 

small   temperature   and   pressure   variations.    Any   random 

change  in  length  in  one  arm  of  the  fiber  interferometer  with 

respect   to   the  other  arm  appears  in  equation  (2.23)   as   a 

change  in  the  angle  (p.    As  (h    changes,  the  amplitudes  of  the 

AC  components  of  equation  (2.23)  varies  due  to   the   factors 

cos  (p   and  sin®.   This  unwanted  change  in  the   amplitude   of 

the   detected   signal   is   referred   to   as   signal   fading. 

Consider  the  case  in  which  the  outputs  of   the   two   arms  of 

the  interferometer  are  initially  in  phase,  CD     is   equal   to 

zero,  and  when  combined  at  the  output  coupler  they  interfere 

constructively.    If   now   there   is   a   small   acoustically 

induced  phase  modulation  e.g.   at  peak  to  peak  amplitude  0.1 

radian j   referring  to  Figure  2.S,  the  photodetector   current 

will  oscillate  about  its  maximum  value  at  twice  the  acoustic 

frequency.   If  the  length  of  one   arm   increases   due   to   a 

thermal    change,   the   average   output   intensity   of    the 
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Figure  2.5    Photodetector  Output  Current  as 

Function  of  Phase  Difference  between 
Interf erometric  Arms 


36 


interferometer  decreases,  but  the  acoustically  induced 
oscillations  will  increase  in  amplitude  and  approach  a  pure 
sinusoidal  o-f  -frequency  equal  to  the  acoustic  -frequency.  In 
a  practical  hydrophone  it  is  essential  to  maintain  the 
interferometer  continuously  at  its  maximum  sensitivity.  The 
sensitivity  to  phase  changes  varies  as  the  slope  o-f  the 
curve  o-f  Figure  2.5.  The  maximum  sensitivity  occurs  -for  odd 
multiples  o-f  7T/2  while  essentially  zero  sensitivity  occurs 
-for  even  multiples  o-f  7T/2.  There-fore  to  maintain 
sensitivity  a  7T/2  phase  bias  is  required.  This  7T/2  bias 
condition  is  known  as  quadrature. 

To  place  the  interferometer  in  quadrature  condition 
for  maximum  sensitivity  and  in  addition  provide  a  means  to 
compensate  for  slow  drifts,  a  phase-locked  compensator 
circuit  is  used.  The  circuitry  required  to  provide  and 
insure  quadrature  operation  in  the  presence  of  low  frequency 
drift  is  shown  schematically  in  Figure  2.6  CRef.  93.  The 
outputs  of  the  two  photodiodes  are  combined  in  a 
differential  amplifier  that  provides  common-mode  rejection 
as  well  as  amplification.  This  is  followed  by  two  stages  of 
integration  that  provide  additional  amplification.  These 
two  integrator — amplifiers  pass  all  signals  from  DC  up  to  the 
highest  frequency  of  interest.  The  output  of  the  two  stage 
integrator-amplifier  is  applied  to  a  phase  shifter  located 
in  the  reference   arm   of   the   interferometer.    The   phase 


37 


unbiased 
photodiodes 


two  stage    integrator 
amplifier 


differentia/ 
amplifier  for 

common  mode 
rejection 


circuit 
reset 

A 


amplifier, 
high  pass 

filter 


to   PIT 


phase 
shifter 


bandpass 
output 


Figure  2.6    Phase-Locked  Loop  Homodyne  Detection 
Circuit 
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shi-fter  consists  of  a  lead  zirconate-lead  titanate 
piezoelectric  (PZT)  cylinder  around  which  the  fiber  in  the 
reference  arm  is  rather  tightly  wound.  The  output  of  the 
integrator — amplifier  is  just  equal  to  the  voltage  required 
to  compensate  for  phase  changes  associated  with  the  low 
frequency  drift  and  that  induced  by  the  acoustic  pressure 
signal  of  interest.  The  effect  then  is  to  produce  a  phase 
shift  in  the  reference  arm  almost  equal  to  that  in  the 
sensing  arm,  causing  the  interferometer  to  remain  very  close 
to  perfect  balance,  i.e.,  to  phase-lock  the  system  CRef. 
103.  If  the  phase  were  exactly  locked  there  would  be  no 
output  signal  from  the  interferometer.  However,  there  must 
be  an  error  signal  or  resulting  phase  modulation  at  the 
photodetectors  in  order  to  have  a  feedback  signal  and  this 
is  effectively  1/G  times  the  modulation  that  would  exist 
without  the  feedback,  where  G  is  the  overall  gain  of  the 
feedback  circuit.  If  the  system  is  initially  at  a  bias 
(operating  or  quiescent)  point  away  from  quadrature,  the 
effective  gain  will  not  be  at  its  maximum  value  and  the 
system  will  tend  to  drift  toward  its  maximum  gain  point  and 
therefore  toward  quadrature.  The  output  signal  of  the 
feedback  system  represents  the  correction  signal  required  to 
keep  the  phase  difference  between  the  interferometer  at  7T/2 
radians.   It  contains  drift  noise  and  signal  information. 

The  signal  out   of   the  compensating  circuit  is  also 
fed  through  a  high-pass  filter  that  has   its   low   frequency 
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limit  set  at  the  lowest  -Frequency  of  interest  but  above  the 
•frequency  o-f  thermal  drift  etc.  There-fore,  the  resulting 
output  is  a  band  o-f  frequencies  corresponding  to  the 
frequency  range  of  interest.  This  constitutes  the  output  o-f 
the  interf erometric  sensor. 

Operational  amplifiers  (OPAMPS)  »rm  used  in  the 
feedback  circuit  and  complementary  metal  oxide  semiconductor 
(CMOS)  components  in  the  reset  circuit.  The  levels  of 
voltage  that  can  be  applied  by  these  circuits  to  the  phase 
shifter  »rm  the  order  of  +  10  volts.  On  the  other  hand,  the 
range  that  the  piezoelectric  phase  shifter  can  accommodate 
is  hundreds  of  volts.  Furthermore,  in  many  cases  the 
amplitude  of  the  phase  drift  resulting  from  temperature  or 
pressure  changes  is  much  larger  than  the  phase  shift  that 
would  be  generated  by  applying  +  10  volts  to  the  phase 
shifter.  Thus,  the  circuitry  keeps  track  of  how  much 
voltage  has  been  applied  to  the  phase  shifter  and  if  the 
limit  of  the  circuit  begins  to  be  reached  it  automatically 
and  rapidly  resets  the  circuit  back  to  its  mid  range  from 
which  point  it  restarts  the  compensation  process.  This  is 
the  purpose  of  the  reset  circuit  indicated  in  Figure  2.6. 
The  phase  change  associated  with  a  large  amplitude  slow 
drift  is  compensated  by  a  number  of  saw  toothed-like  small 
amplitude  phase  changes.  Care  must  be  taken  in  the 
alignment  of  the  compensation  circuit  to  minimize  noise, 
introduced  during  the  reset  process. 
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"I-   EXPERIMENTAL  APPARATUS 
A.   INTERFEROMETER  DESIGN  AND  CONSTRUCTION 

Figure  3.1  shows  a  schematic  diagram  of  the  experimental 
interf erometric  system  utilized  in  this  study.  The 
interferometer  is  modular  in  construction.  As  each  module 
was  completed  proper  operation  was  verified.  A  brief 
description  of  the  overall  construction  of  the 
interferometer  is  first  presented;  this  is  followed  by  a 
detailed  description  of  each  individual  part. 

The   system   is  comprised  of  an   all -fiber   Mach-Zehnder 

interferometer.   A  Laser  Diode   Labs   solid   state  laser  was 

1 
pigtailed    to  high  numerical  aperture  single-mode  (830   nm) 

Corning  fiber.   The  laser,  with   the   pigtailed   fiber,   was 

spliced  (see  section   on   fiber  preparation  and  splicing)  to 

approximately  half  a  meter   of   ITT   single-mode  fiber.   The 

free  end  of  the  ITT  fiber  was   then   spliced  to  an  ITT  2X2 

all -fiber  3dB  coupler.   Initial  attempts  were  made  to  splice 

the  pigtailed  Corning  fiber   directly   to   the  coupler  fiber 

but   with   limited   success.    Differences   between   melting 

temperatures -of  the  Corning  fiber  and  coupler  fiber  resulted 


1  Pig tailing  of  fiber  to   laser   accomplished   by   laser 
manufacturer 
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I.  Laser  diode  5. Sensing  coil 

2.3dB  Coupler  6.  Photodetector 

3  Polarization  controller 
4.PZT  Phase  shifter       X-Indicates  splice 


Figure  3.1    Schematic  of  Experimental 
Inter-ferometric  System. 
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in  low  strength,  high  loss  splices.  For  this  reason  the  ITT 
T-1601  -fiber  had  to  be  used  as  an  intermediate  link  between 
the  Corning  -fiber  and  ITT  coupler  -fiber. 

The  couplers  themselves  are  all  -fiber  counterparts  o-f 
bulk  optical  beam  splitters  and  equally  split  the  light  at 
the  input  end  of  the  interferometer.  Half  the  light  is 
coupled  to  one  arm  and  the  remaining  light  into  the  other 
arm. 

Two  single-mode  optical  fibers  with  approximately  equal 
length  (estimated  path  length  difference  of  approximately  2 
cm)  form  two  arms  of  the  interferometer.  Both  were 
constructed  identically;  therefore  only  the  construction  of 
one  will  be  discussed.  Eight  meters,  from  a  ten  meter 
length  of  fiber,  was  wrapped  in  a  single  layer  onto  a 
piezoelectric  (PZT)  cylinder.  At  each  end  of  the  wrapping 
approximately  a  meter  of  lead  was  left.  One  lead  was  then 
spliced  to  fifteen  meters  of  fiber  of  which  approximately 
ten  meters  was  coiled  into  a  4  cm  diameter  coil.  The  other 
lead  coming  from  the  fiber  wrapped  around  the  PZT  is  fed 
through  mechanical  polarization  controllers  CRef.  113.  One 
end  is  then  spliced  to  one  of  the  output  leads  of  the  input 
coupler  and  the  other  end  is  spliced  to  one  of  the  input 
leads  of  the  output  coupler  to  complete  construction  of  an 
arm  of  the  interferometer.  The  output  coupler  recombines 
the  light  from  the  two  arms.  Basically  the  interferometer 
is  used   to  convert  phase  modulation  to  intensity  modulation 
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which  occurs  when  light  -from  the  two  arms  is  combined.  The 
two  outputs  of  the  interferometer  are  detected  by  two 
phot odi odes. 

Figure  3.2  is  a  photograph  o-f  the  completed 
interferometer.  The  design  of  the  system  allows  for  a 
single  coil  to  be  placed  in  water  to  act  as  a  pressure 
hydrophone  or  for  both  coils  to  be  placed  in  water  for 
operation  as  a  pressure  gradient  hydrophone. 

B.   INTERFEROMETER  COMPONENTS 

1-   Laser  Source 

The  optical  source  used  in  this  experimental 
interf erometric  system  was  a  single  frequency  infrared  laser 
diode  optimized  to  emit  coherent  light  at  a  wavelength  of 
830  nanometers.  Two  solid  state  laser  diodes  from  Laser 
Diode  Labs  were  used.  The  first  was  a  CW-234  laser, 
pigtailed  by  the  manufacturer  with  a  short  section  of  high 
numerical  aperture  single  mode  (830  nm)  Corning  fiber.  The 
second  laser,  a  CWX-233  laser  also,  pigtailed  with  Corning 
fiber,  replaced  the  first  laser  after  an  inadvertent 
increase  in  current  destroyed  the  CW-234  laser.  Power  for 
the  laser  was  provided  by  a  12  volt  18  amp-hr  Gel  cell 
battery  through  a  voltage  regulating  circuit.  This  circuit 
maintained  constant  current  to  the  laser  diode.  Figure  3.3 
is  a  schematic  of  the  power  supply.  Before  using  either 
laser  in  the  interferometer  the  end  of  the  Corning  fiber  was 
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Figure  3.2    Photograph  of  Experimental 
Interf erometric  System 
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Figure  3.3    Schematic  of  Laser  Diode  Power  Supply 
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cleaved  and   set   up   to  enable  a  Coherent  212  optical  power 

meter  to  sense  the  laser  diode  power  output.    As  current  to 

the  laser  diode  was  increased  from   0   to   approximately   53 

mi  11  lamps  -For   the   CW-234   laser,  and  0  to  39  mi  11  lamps  for 

the  CWX-233  laser,  the  power  of  light  emitted  from  the  laser 

and  coupled  into  the  fiber  was   read   on   the   power   meter. 

Figures  3.4  and  3.5  Are   graphs   of   laser   power   vs   source 

current   for   the  CW-234  and  CWX-233   lasers,   respectively. 

Extrapolating  the  linear  region  back  to  zero  power  yields   a 

threshold  current  of  48.5  mi  11  lamps  for  the  CW-234  laser  and 

31.3  milliamps  for  the  CWX-233  laser.   The   current   to   the 

laser  diode  was  maintained  at  53  milliamps   for   the  CW-234 

laser   and    35   milliamps   for   the   CWX-233   laser   during 

experimental  testing  and  data  acquisition. 

2-   Ei^?C  Sp_eci.f  i9.ati_gn 

The  optical  fiber   used   was   ITT  single-mode  fiber, 
2 
type   T-1601   optimized  at  0.83  microns.   Specifications   of 

fiber  are  as  fallows: 

Fiber  I dents  81092U-17Bic 

Preform  No.:  EMT-21972B 

Numerical  Aperture:  0.12 

Core  Diameter:  4.57  microns 


2  Optical  fiber  obtained  from  International  Telephone 
and  Telegraph,  Electro  Optics  Products  Division,  Roanoake, 
Va  24019 
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Outer  Cladding  Diameter:  75  microns 

Primary  Sheath:  6E  615  Silicone 

Secondary  Sheath:  polyester  Hytrel 

Total  Diameter:  16  mil 

Attenuation:  2.07  dB/Km  at  0.83  microns 

Figure  3.6  is  a  sketch,  provided  by  the   manufacturer  of  the 
end  o-f  the  fiber,  with  typical  dimensions,  approximate  index 
o-f  re-fraction  profile,  and  typical  spectral  attenuations. 
3.   3  dB  Coupler 

The  2  X   2   single-mode  couplers  were  also  purchased 
from  ITT.  The  following  specifications  are  provided: 

Serial  Nos. :  JM-SM-108  &  JM-SM-109 

Fiber  No.:  EMC  41556C/830427-401a 

Fabrication  Date:  6/3/83 

Excess  Loss:  0. 1  dB 

Uniformity:  0.1  dB 

Operating  Wavelength:  0.83  microns 

The  purpose  of  a  3  dB  coupler  is  to  split  the  light  equally 
into  the  two  arms  of  the  interferometer  or  to  recombine  the 
light  causing  the  light  from  the  two  arms  to  interfere  in 
the  output  fibers  and  on  the  face  of  a  photodiode. 
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DIMENSIONS   SHOWN   ARE   NOMINAL  VALUES 
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TYPICAL  SPECTRAL  ATTENUATION  -  SINGLE  MODE  OPTICAL  FIBER 
Figure    3.6  ITT    Single    Mode    Fiber    T-1601 
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4"   ElSSgelectric  Phase  Shifter 

The  phase  shifter  consists   of  a  lead  zirconate-lead 
3 
titanate  (PZT)  cylinder    around   which  the  fiber  is  tightly 

wrapped.   The  cylinders  used   were   type   Channel  5500.  Each 

cylinder   is   3.8   cm   (1.5   in.)   long  by  3.8  cm   in   outer 

diameter  with  a  wall  thickness  of  0.32  cm   (0.225  in.).   Two 

of    these   cylinders   were    incorporated     into    the 

interferometer,  one  in  each  arm.   The  PZT  cylinders   have   a 

relatively  small  dynamic  range  but  by  wrapping  many  turns  of 

fiber  on  the  PZT  it   was   possible  to  increase  the  magnitude 

of  the  phase  shift.   Chapter   II,   Section   C   presents   the 

theory  of  operation  of  the  PZT  as  a  phase  shifter. 

5.   PgL»Ci?*tign  Controller 

A  polarization  controller,  as  described  by  Lefevre 
CRef.  113,  was  used  in  this  experiment.  A  very  simple 
configuration  is  employed  to  control  the  polarization  of 
single  mode  fiber.  This  device  is  equivalent  to  fractional 
wave  plates  of  classical  optics.  The  controller  uses  the 
stress  birefringence  induced  by  bending  the  fiber. 

°-   Eb-Qt-Qd^tectgr 

The  photodetector  used  to   detect  the  optical  output 

4 
of   the   fiber   interferometer   is   a   Clairex   type   CLD-41 


3  Piezoelectric   cylinders    obtained    from    Channel 
Industries,  Inc.,  Santa  Barbara,  California  93105 

4  photodiodes  obtained  from  Clairex  Electronics,   Mount 
Vernon,  New  York,  New  York  10550 
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photodiode.  They  are  all  silicon  PN  planar  diodes  with  high 
linearity,  low  dark  current,  and  fast  response.  The 
•following  are  the  electrical  characteristics  as  given  by  the 
manufacture: 

Active  Area:  1.3x1.3  mm 

Short  Circuit  Current:    min  6  to  max  12  microamps 
Qpen  Circuit  Voltage:     .40  volts  typical 
Dark  Current:  1  microamp 

Junction  Capacitance:     200  picofarads 
Rise  or  Fall  Time:        5  microseconds 
Temperature  Coefficient:  +.2%/degree  C  typical 
Peak     Spectral     Response:   0.91  microns 
Figure  3.7  is  a  plot  of   photodiode   DC   voltage   vs   system 
laser  optical  power  as  measured  in  the   laboratory   using   a 
Hewlett-Packard  HP-3438A  digital  multimeter  to   measure   the 
output  voltage. 

C.   FIBER  PREPARATION  AND  SPLICING 

Before  splicing  between  fibers  could  be  accomplished, 
each  end  of  the  fibers  had  to  be  prepared  by  stripping  and 
cleaving.  Fiber  stripping  was  accomplished  in  two  steps. 
First,  the  outer  Hytrel  plastic  sheath  was  removed  by  use  of 
a  razor  blade.  The  fiber  was  placed  on  a  piece  of  black 
construction  paper  and  then  a  razor  placed  at  a  very  small 
angle  to  it.  With  a  slight  pressure  applied  to  the  razor, 
the  fiber  was  drawn  out  from   underneath   the   razor.    This 
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Mill  carve  the  Hytrel  jacket  neatly  -from  the  -fiber.  By 
rotating  the  fiber  and  repeating  this  process  several  times, 
most  of  the  outer  Hytrel  and  inner  silicone  sheathing 
material  will  be  removed. 

The  second  step,  that  of  removing  any  residue  of  inner 
or  outer  sheath  that  was  not  removed  by  the  razor  blade,  is 
best  done  chemically  with  a  solution  of  chromic  acid.  This 
is  done  by  immersing  the  section  of  fiber  that  has  been 
partially  stripped  by  the  razor  blade  in  the  acid  for 
approximately  two  minutes.  Any  jacket  material  that  remains 
will  turn  to  a  jelly-like  consistency  which  is  easily 
removed  by  wiping  with  a  soft  tissue.  After  rinsing  in 
distilled  water,  the  fiber  is  pulled  between  a  tissue 
saturated  with  methyl  alcohol.  If  the  fiber  is  completely 
clean  of  the  jacket  a  squeaking  sound  will  be  heard.  At 
this  point  the  fiber  is  bared  to  the  silica  cladding  and 
must  be  handled  with  caution.  An  unscratched  fiber  has  very 
high  tensile  strength,  but  will  break  easily  if  scratched 
and  bent. 

In  order  to  insert  and  extract  light  efficiently  to 
obtain  a  low  loss  splice,  the  ends  must  be  flat  and 
perpendicular  to  the  longitudinal  axis  of  the  fiber.  The 
process  utilized  to  achieve  this  was  cleaving.  To  cleave 
the  fiber,  the  ends  are  first  stripped  back  about  two 
centimeters  longer  than  needed.  After  preparation  for 
cleaving,   the   fiber   was   held   with   slight   tension   and 
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carefully   touched   to    the    cleaving   blade   at   a   point 

5 
approximately  2  cm  from  the   end   of  the  fiber.   The  cleaved 

^rrnrn   was   then   placed  on  a  smooth,  round  mandrel  that  has  a 

radius   of   approximately  four  centimeters.   The   fiber   was 

held  with  thumb  of   one   hand.    Maintaining  modest  tension, 

the  fiber  was  pulled  around  mandrel  with  other   hand.    When 

scored  »rmai      reaches  mandrel,  a  clean   break   will   usually 

result.   Experience  showed  that  if   the   fiber   broke   while 

being  scored  it  would  have  a  bad  end  face. 

To  check  the  cleave,  the  fiber  is  inspected  under  a 
microscope.  Angular  breaks  other  than  90  degrees  will  not 
make  a  good  low  loss  splice.  In  addition,  chips  and  barbs 
at  the  cleaved  ends  will  not  allow  satisfactory  splices. 
Figure  3.8  shows  a  good  cleave  along  with  improper  cleaves. 

Splicing  of  fiber  was  accomplished  using  an  electrical 
arc  welder.  The  type  of  welder  used  was  a  Power  Technology, 
Inc.  portable  fiber  splicer  model  PFS-240.  Basic  splicer 
operations  of  pref using  and  fusing  for  monomode  fiber  were 
done  in  accordance  with  the  operating  manual.  Thumbwheel 
settings  which  controlled  splicer  arc  amperage  and  the  time 
variation  of  the  arc  current  usually  used  »rB   as  follows: 

Ramp  Time     Arc  Time     Arc  Current 
0.1  sec      0.6  sec      16  mA 


5  Sapphire   cleaver   obtained  from  General  Fiber  Optics, 
Caldwell,  New  Jersey  07006 
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Final  alignment  o-f  cores  -for  maximum  light  transmission 
was  used  to  produce  a  low-loss  splice.  The 
micro-positioners  on  the  splicer  were  adjusted  to  bring  the 
prepared  ends  o-f  the  -fiber  into  close  alignment.  Light  from 
the  laser  was  injected  into  the  -free  end  o-f  one  o-f  the 
fibers,  transferred  across  the  proposed  splice  to  the  second 
fiber,  and  the  output  from  the  latter  was  projected  onto  the 
face  of  the  power  meter's  photodiode.  The  mi cropositi oners 
were  adjusted  until  the  power  meter  indicated  that  the 
transmission  across  the  splice  was  maximum.  At  that  point 
the  fiber  ends  were  brought  together,  heated  by  the  arc 
welder  and  fused.  As  the  joint  cools,  surface  tension  will 
pull  the  fibers  into  final  alignment.  Reports  in  the 
literature  indicate  that  losses  obtained  by  fusion  splicing 
may  be  as  low  as  0.1 's  dB.  However,  in  this  study  many 
hours  of  practice  were  required  to  obtain  good  (i.e.  less 
than  3dB  loss)  splices.  Nine  splices  are  required  in  the 
configuration  of  the  experimental  system  described  above. 
The  average  loss  from  each  splice  was  approximately  ldB  and 
the  total  loss  of  the  constructed  interferometer  was 
approximately  8.4dB. 

Prior  to  projecting  light  from  the  end  of  the  fiber  onto 
the  surface  of  a  photodiode,  the  end  of  the  fiber  was 
prepared  in  the  same  manner  as  for  splicing  with  the 
exception  that  the  length  of  fiber  that  had  its  sheath 
removed  was   painted  black.   Painting  was  used  to  accomplish 
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mode  stripping.  Cladding  light  was  stripped  prior  to 
reaching  the  end  of  the  -fiber,  thus  only  core  light  was 
detected  by  the  photodiode.  In  splicing  each  additional 
fiber  element  into  the  system  it  was  necessary  to  mode  strip 
its  output  end  just  before  the  energy  meter  to  ensure  that 
only  core  light  is  detected,  i.e.  to  absorb  unwanted  light 
that  may  reach  the  output  end  via  the  cladding. 

D.   INSTRUMENTATION  AND  DATA  ACQUISITION 

Basic  block  diagrams  of  the  experimental  systems  used  to 
obtain  data  on  interferometer  performance  will  be  shown  in 
Chapter  IV  as  the  experimental  procedures  are  discussed. 
Computer  data  acquisition  was  used  for  portions  of  the  data 
taking.  The  Hewlett-Packard  HP-85F  computer  coordinates  the 
peripherals  and  records  and  displays  the  data.  The  printer 
and  plotter  are  non-essential  to  the  system  and  were  used  to 
provide  larger  format  graphics  and  printouts  than  are 
available  with  the  built-in  thermal  printer  provided  with 
the  HP-85F.  The  fallowing  is  a  brief  description  of  each 
individual  piece  of  equipment. 

The  HP-85F  is  an  eight  bit  microprocessor  that 
utilizes  BASIC  computer  language.  The  computer  has  as 
standard  16K  bytes  of  read/write  memory.  16K  bytes  of 
additional  memory  were  added  to  the  system  giving  a  total  of 
32K  bytes.   The  HP-85F  has  a  127  millimeter   diagonal   black 

59 


and  white  electromagnetic  CRT.  A  32  character  per  line 
thermal  printer /pi otter  is  part  of  the  unit.  Programs  or 
data  may  be  stored  on  and  read  from  magnetic  tape 
cartridges.  To  interface  with  peripheral  equipment,  an  I/O 
ROM  and  an  interface  card  were  added  to  provide  HP-IB  (IEEE 
standard  488-1975)  instrumentation  capabilities. 

2.   Arbi.t rary_  Waveform  Generator  HP  3314A 

The  Hewlett-Packard  Model  3314A  arbitrary  waveform 
generator  can  produce  three  different  pre-programmed 
waveforms)  sine,  square,  and  triangular.  This  instrument 
has  a  frequency  range  of  1  microHertz  to  20  megaHertz. 
Output  amplitude  is  0  millivolts  to  10  volts  peak  to  peak 
into  a  fifty  ohm  load.  This  model  is  fully  programable 
through  a  HP- IB  connection.  For  this  experiment  sine  waves 
of  various  frequencies  and  amplitudes  were  usually 
generated.  During  alignment  of  the  locking  circuit, 
however,  a  small  amplitude  100  Hertz  square  wave  was  used. 

3-   §e?ctrums  Analyzer  HP  3582 A 

The  Hewlett-Packard  Model  3582A  is  a  dual  channel 
spectrum  analyzer.  This  instrument  has  a  frequency  range  of 
0.02  Hertz  to  25,600  Hertz.  The  analyzer  has  a  11.9  by  9.6 
cm  CRT  that  can  display  two  simultaneous  information  traces, 
plus  four  lines  of  alphanumeric  data  giving  measurement 
configuration  and  results.  Frequency  spans  from  1  Hertz  to 
25,000  Hertz  full  scale  allow  flexibility  in  selecting  the 
portion  of  the  spectrum  to  be  analyzed.   Spans   from  5  Hertz 
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to  25,000  Hertz  can  be  positioned  anywhere  within  the 
•Frequency  range  o-f  the  instrument  to  provide  excellent 
-frequency  resolution.  The  instrument's  "front-end" 
sensitivity  ranges  can  measure  and  analyze  -from  1  micro-volt 
to  31.6  volts  and  has  a  dynamic  range  of  70  dB. 
4.   Qsci.l.l.gscgpe  lektrgnix  2215 

The  Tektronix  Model  2215  oscilloscope  is  a  dual 
channel,  60  MHz  instrument.  Its  vertical  system  provides 
calibrated  deflection  factors  from  2  millivolts  to  10  volts 
per  division.  Trigger  circuits  enable  stable  triggering 
over  the  full  bandwidth  of  the  vertical  system.  The 
horizontal  system  provides  calibrated  sweep  speeds  from  0.5 
seconds  per  division  to  50  nanoseconds  per  division. 

5-  Digital  Multimeter  HP  3438  A 

This  multimeter  was  utilized  to  obtain  rms  voltage 
readings  of  an  LC-10  hydrophone  output  and  transmit  this  via 
the  HP- IB  to  the  HP-85  computer. 

6-  Standard  Hy.drgp.hgne  LC-10 

6 
An  LC-10  hydrophone  ,  Serial   No.   2341   was  used  as 

the  standard  hydrophone  for   sensitivity  calculations  of  the 

single  coils  and  of  the  dual  (gradient)  coil  setup.   Average 

free   field   voltage   sensitivity   for   this   hydrophone   is 

specified  by  the  manufacture  to  be  -209.2  dB  re  1  Volt//!  Pa. 


6  Hydrophone   manufacture   is  Celesco  Ind. ,  Canoga  Park, 
California,  91304 
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This  hydrophone  and  an  USRD  type  F27  reversible  transducer 
were  used  in  another  experiment  CRef.  123  to  obtain  a 
reciprocity  calibration  of  an  USRD  type  F33  transducer.  The 
sensitivity  -for  the  LC-10  hydrophone  obtained  from  the 
experiment  was  within  2  dB  of  manufacturer's  specifications 
between  5  an  30  KHz. 

7.   Acoustic  Calibrator 

An  acoustic  system  was  constructed  for  calibrating 
hydrophones  of  various  types  by  the  comparison  method. 
Figure  3.9  shows  the  calibrator  in  the  lab,  and  Figure  3.10 
shows  the  details  of  construction.  The  calibration  section 
consists  of  an  open-ended  column  of  water  excited  by  a 
moving-coil  driver.  A  USRD  type  J 11  transducer  is  used  as 
the  driver.  The  calibration  chamber  consists  of  a  6  inch 
inner  diameter,  1/4  inch  thick  PVC  1120  tube  (ASTM  1785), 
mounted  on  a  rim  around  the  face  of  the  driver. 
Compensation  for  the  hydrostatic  pressure  of  the  water  on 
the  face  of  the  driven  piston  is  provided  by  the  airtight 
enclosure  mounted  under  the  driver  assembly.  A  manometer  is 
used  to  measure  the  equalizing  pressure.  A  valve  is  opened 
and  air  is  pumped  into  the  equalizing  chamber  until  the 
internal  air  pressure  is  equal  to  the  water  pressure  on  the 
upper  face  of  the  driver.  This  system  is  similar  in  design 
to  the  USRD  type  640  shipboard  calibrator  used  by  the  U.S. 
Navy  CRef.  13D  ,  based  on  article  by  Schloss  and  Strasberg 
CRef.  143. 
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Figure  3.9    Photograph  of  Acoustic  Calibrator 
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Figure  3.10    Sectional  View  of  Acoustic  Calibrator, 
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The  compliance  of  the  PVC  tube  reduced  the  speed  o-f 
sound  o-F  the  water  in  the  calibrator  to  approximately  1/5  o-F 
the  normal  speed  o-f  sound  in  water.  This  was  useful  in  the 
experiments  since  it  reduced  the  wavelength  A  by  a  factor 
of  5  at  any  given  frequency  and  resulted  in  a  larger 
gradient.  Figures  3.11,  3.12,  and  3.13  show  the  standing 
wave  acoustic  field  at  three  resonant  frequencies  of  the 
calibrator.  The  figures  are  graphical  plots  of  LC-10 
hydrophone  voltage  as  a  function  of  position  in  the 
calibrator.  By  measuring  the  distance  of  a  half  wavelength 
at  the  three  resonant  frequencies  the  speed  of  sound  in  the 
calibrator  can  be  calculated.  From  the  data  shown  in  Table 
I  an  average  value  of  358.3  meters  per  second  with  a 
standard  deviation  of  8.6  m/sec  for  the  speed  of  sound  in 
the  calibrator  is  calculated.  The  data  was  gathered  over  a 
period  of  one  month. 

TABLE  I 


Sound  Speed  in  Calibrator 


Frequency  (Hz) 

390 
390 
391 
595 
593 
593 
593 
593 
752 


Half 
Wavelength 
44.5 
44 

45.5 
30 
31 

30.5 
30.5 
30.5 
24.5 


(cm) 


Speed  of 
Sound  (m/sec) 
347.1 
343.2 
355.8 
357.0 
367.7 
361.7 
361.7 
361.7 
368.5 


average  >>>  358.3  +8.6 
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Figure  3.12    LC-10  Voltage  in  Calibrator  at  593  Hz 
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8-   Locking  QiC£ui_t 

A  phase  lock  loop  locking  circuit  similar  to  the 
circuit  presented  by  Fritsch  and  Adamovsky  CRef.  93  and 
discussed  in  Chapter  II  was  used  to  verify  operation  of  the 
interferometer  as  a  single  and  dual  (gradient)  hydrophone 
system.   Data  will  be  presented  in  the  next  chapter. 
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IV.   EXPERIMENTAL  PROCEDURE  AND  RESULT;? 
A.   INTERFEROMETER  CHARACTERISTICS 

A-fter  construction  of  the  inter -f erometer  was  completed, 
testing  was  conducted  to  verify  proper  operation  of  the 
interferometer  as  a  sensor.  The  initial  tests  consisted  of 
detailed  measurements  taken  to  determine  the  amplitude  of 
the  optical  phase  shift  as  a  function  of  the  drive  voltage 
applied  to  the  piezoelectric  cylinder  in  one  of  the  arms  of 
the  interferometer.  Figure  4.1  shows  a  block  diagram  of  the 
instrumentation  required  for  operating  the  system.  A 
variable  amplitude  sine  wave  of  frequency  f  from  the 
function  generator  was  applied  to  the  PZT  cylinder.  A 
Tektronix  2215  oscilloscope  and  a  HP  3582A  spectrum  analyzer 
were  used  to  detect  the  output  of  one  of  the  photodiodes. 

The  behavior  of  a  Mach-Zehnder  fiber  interferometer  was 
treated  in  detail  in  Chapter  II.  Recall  that  the  AC  portion 
of  the  photodiode  current  is  given  by  the  following 
equations 


<fi   *  2  £ 


CO 


i  (t)  -  i  cos<t)  *  2  )    J       (k  x)cosC2n(  (jj  t)3 

0    -  L  n^  2n   o  s 

-i  sin<j5)*  2^J     <k  x)sinC  <2n+l)  (  Gl)  t>  3 
0       L  pj-Q   2n+l   o  s 


where  J   is  the  Bessel  function  of  order  n. 
n 
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Figure  4 . 1 


Block  Diagram  of  Instrumentation  Setup 
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From  equation  (4.1)  it  is  clear  that  as   the  phase  shift 

k  x  is  increased   the   photodiode   output  becomes  non-linear 

0 
and  harmonics  o-f  the  drive  signal  begin  to   appear.   This  is 

illustrated    in   Figure   4.2   where   oscillograms   o-f    the 

photodiode   output  signal  are  shown  -for  successively   higher 

drive   signals   applied  to  the  piezoelectric   cylinder.    In 

each  oscillogram  the  drive  signal   is   shown   in   the  upper 

trace  and   the   corresponding   photodiode   output   signal  is 

shown  in  the  lower  one. 

The  magnitude  o-f  the  optical  phase  shift  as   a   function 

of   drive   voltage   may    be   determined   quantitatively   by 

examining  the  frequency   spectrum   of   the  photodiode  output 

signal.   From   equation   (4.1)   it   may   be   seen   that   the 

th 
amplitudes  of  the  various   n     order   harmonics  should  have 

nulls   for   phase   shifts   k  x   equal   to  the  zeroes  of   the 

0 
th 
corresponding  n    order  Bessel   functions.   This  behavior  is 

illustrated  in  Figure  4.3   where   the   photodiode  voltage  is 

plotted  for  the  fundamental  and  first  two  harmonics  vs  drive 

voltage   applied   to  the  piezoelectric  cylinder.    The   data 

were   obtained  with  the  computer  data  acquisition   equipment 

discussed  in  chapter  III.  A  computer  program  for  the   HP-85F 

computer,   listed   in   Appendix   A,   was   written   in   BASIC 

computer  language  to   control   the   peripheral  equipment.   A 

fundamental  frequency,  lower  and   upper  limits  of  amplitude, 

and   an   incremental   value    are   selected.    The   function 

generator  and  spectrum  analyzer  are  preset  and  configured  to 
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A  / 


0.5  volts  p-p   drive 
Chi:  0.2vott/div 


A    A    A    A 

y  V  v  SJ  \ 


Ch  2:   2  m  volt/div 


2  volts   p-p  drive 
Chi:    I  volt/div 


A   A  A   A  A   A  A 
I  /  V  \  l\J\  /  V 


C7?2:  2  m  volt/div 


3  volts  p-p    drive 
Chi:   I  volt/div 


\ 


Ch2:  2mvoit/div 


Figure  4 . 2 


Oscillograms  of  PZT  Drive  Signal  and 
Photodiode  Output  Signal 
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the  prescribed  settings.  At  the  loner  voltage  amplitude 
limit  the  spectrum  analyzer  seeks  -for  the  maximum  amplitude 
o-f  the  fundamental  and  -first  two  harmonics  attained  during  a 
selected  period  of  time.  The  amplitude  of  these  signals 
vary  in  time  as  the  angle  (7)  varies  randomly.  The  sample 
time  selected  on  the  spectrum  analyzer,  of  the  order  of  16 
seconds,  was  sufficient  that  the  actual  maximum  values  were 
attained  and  recorded.  The  HP-85F  computer  directs  the 
spectrum  analyzer  to  output  the  maximum  amplitude  values  of 
the  fundamental  and  harmonics.  These  values  are  stored  by 
the  computer.  The  computer  then  directs  the  function 
generator  to  increase  the  drive  voltage  by  the  specified 
incremental  value.  The  spectrum  analyzer  again  seeks  for 
the  maximum  amplitudes  and  the  process  is  repeated.  This  is 
done  for  each  incremental  value  until  the  upper  limit  is 
reached.  The  data  stored  in  the  computer  memory  is  sent  to 
the  printer  and  also  transferred  to  magnetic  tape  for 
permanent  storage  to  enable  analysis  at  a  later  time.  A 
graph  of  the  data  is  plotted  on  the  computer  printer. 
Sample  of  data  obtained  is  seen  in  Appendix  B.  From  this 
data  plots  of  the  type  shown  in  Figure  4.3  could  be 
generated. 

The  computer  program  was  run  initially  with  the  CW-234 
laser  source.  Three  runs  were  made,  two  with  the  fiber 
wrapped  PZT  cylinder  of  one  arm  acting  as  the  modulator  and 
one  with  the   other   fiber   wrapped   PZT   of   the  other  arm. 
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After  replacement  of  the  source  with  the  CWX-233   laser   was 

accomplished  two  more   runs   were   made,   one   on   each   PZT 

cylinder.   From  the  data,  the  amplitude  of  the  drive   signal 

that  causes  the  fundamental  to  be   at   its   first   null   was 

determined  to  be  1 . 65V (peak) +0. 04V (peak) .  Thus,  referring  to 

equation  (4.1),  it  may  be  concluded  that  the  drive  condition 

was  such  that  the  first  order  Bessel  function   J    was  at  a 

1 
zero.   The  first  of  the  fundamental 's   zeroes   occurs   at   a 

phase  shift  of  k  x  =  3.83   radians.    Therefore   for   the   8 

0 
meters  of  fiber  wrapped  on  the  PZT  cylinder  a  value  of  0.29 

rad/volt (peak) /meter    results.    By   observing   the    drive 

voltages   required   to    produce    the   first   zero   of   the 

fundamental   as   the  frequency  was  varied   manually   it   was 

possible  to  determine  that  the  phase  shift  per  volt  remained 

constant,  within   the   precision   of   the   measurements,  for 

frequencies  from  80  hertz  to  1500  hertz. 

Agreement  of  experimental  data  presented  in  Appendix  B 
and  illustrated  in  Figure  4.3  to  Bessel  function  theory  is 
shown  in  Table  II.  The  ratio  of  experimental  maxima  and 
minima  of  the  same  harmonic  are  compared  to  the  theoretical 
ratio  of  maxima  and  zeros  of  Bessel  functions  of  the  same 
order.  Theoretical  data  as  tabulated  by  Abramowitz  and 
Stegun  CRef.  153  was  used. 

The  last  test  conducted  was  to  verify  operation  of  the 
system  with  the  locking  circuit.  After  aligning  the 
locking,  a  small  amplitude  sine   wave   was  applied  to  one  of 
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the  PZT  cylinders  simulating  an  acoustic  signal.  The 
locking  circuit  produced  a  signal,  similar  to  this  input 
signal,  that  was  fed  back  to  the  PZT  cylinder  in  the  other 
arm  to  act  as  a  compensation  signal,  thus  maintaining  the 
two  arms  in  quadrature.  This  output  signal  -from  the  locking 
circuit  was  also  sent  to  the  oscilloscope  to  show  that  the 
output  signal  did  -Follow  the  input  signal. 

B.   INDIVIDUAL  SENSOR  SENSITIVITY 

The  individual  acoustic  sensitivity  o-f  the  sensor  coils 
in  each  arm  of  the  interferometer  was  calculated  from 
measurements  obtained  in  the  acoustic  calibrator.  A  single 
coil  was  mounted  on  a  stand  with  an  LC-10  hydrophone  placed 
in  the  center  of  the  coil  as  shown  in  Figure  4.4.  The 
acoustic  calibrator  was  driven  at  a  frequency  that  generated 
a  standing  wave  acoustic  resonance.  Chapter  III  section  D7 
discusses  the  standing  waves  generated  in  the  calibrator. 
The  coil /LC-10  hydrophone  setup  was  lowered  to  the  first 
maximum  point  of  the  standing  wave.  Output  of  one  of  the 
photodiodes  of  the  interferometer  was  sent  to  the  Tektronix 
2215  oscilloscope  and  the  HP-3582A  spectrum  analyzer. 
Output  of  the  LC-10  hydrophone  was  sent  to  the  HP-3438A 
Multimeter  set  to  read  AC  voltage.  A  schematic  of  the 
instrumentation  setup  is  shown  in  Figure  4.5.  The  spectrum 
analyzer  was  set  to  monitor  the  fundamental  frequency  and 
the  first  two  harmonics.   Drive  voltage  was  increased  to  the 
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Figure  4.4     Photograph  of  LC-10  Hydrophone  and 
Single  Fiber  Coil  Setup 
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Figure  4.5     Block  Diagram  of  Calibration  Circuit 
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acoustic   calibrator   until  the  voltage  of   the   -fundamental 

frequency,  as  seen  on  the   spectrum   analyzer,   reached   its 

first   minimum   value.     The    voltage   of   the   LC-10   was 

recorded.   If  there  was  more  than  one  pressure   maximum   the 

sensor  arrangement  was  moved   to   the   next  pressure  maximum 

and  the  LC-10  voltage  at  the  fundamental 's  first  zero   point 

was  recorded.   Data  were  recorded  for  all   of   the   pressure 

maxima   of   the   standing   wave.    Data   were  obtained   both 

manually  and  by  computer  control  as  discussed  in   the   above 

section. 

The  first  minimum  of  the  fundamental  occurs  at   a   phase 

shift  of  3.83  radians,  therefore  the  sensitivity  of  the  coil 

M    in   optical   radians   per  micropascal  is   calculated   as 

C 
follows: 


M  -  M  *3.83rads/V 
C    H  H 


(4.2) 


where  M   is  the   sensitivity   of   the   LC-10   hydrophone   in 

H 
vol ts/Zi Pa   obtained  from  manufacture  specification   sheets. 

V   is  the  LC-10  hydrophone   voltage   where  the  first  zero  of 

H 
the  fundamental  occurred.    Table   III  is  a  summarization  of 

all  single   coil   sensitivity   measurements   taken   at  three 

different   calibrator    resonances    and    different   depths 

associated  with  pressure  maxima. 
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Table  III 
Individual  Coil  Sensitivity 


Coil  #  1 


Coil  #  2 


Freq 'Depth  of 

(Hz) IHyd  from 

! Air-Water 

! Boundary 

i   (cm) 
■ 

iNos. 
!  of 

!  Runs 

1 

: : 

! Avg  ! Std  i  S  Freq 
S  M   Idev  S  S  (Hz) 

>U r ad  i/ir ad!  ! 
\U   Pa !  U  Pa  1  S     I 

1  -3  ■  -3  ! ! 

110   110   II 

Depth  of 
Hyd  from 
Air-Water 
Boundary 
(cm) 

Nos. !Avg  iStd 
of  !  M   idev 

Runs  l/irad  \/J.r*d 
\/J  Pai^Pa 
•  -.3  '  -*? 

!10   !10^ 

392  S    23 

385  !    60 

385  !    80 
■ 

:  l 
:  l 
:  9 

S7.07!   -  ! 

112.41   -  1 

!9. 1512. 31! 
■     i     i 

390 
389 
385 

23 
65 
78 

l 
l 

9 

!5.65i   - 
!9.99!   - 
!10.3!2.4 

■ 

544*S    16 

593  !    48 

544  S    47 
■ 

:  6 

!   7 

:  10 

i     i     i 
•  4 . 95 1 . 563 1 
18.6611. 14 i 
! 9. 93 5  2. 97! 

1        !        ■ 

593 
593 

16 

1    48 

4 
3 

■  i 
i5.19i.314 

! 1 1 . 5  i . 529 

•   _  ■   _ 

■  i 

■  i 

■ 
695* t    12 
695  !    37 
695  !    65.5 
695  !    91.5 

:  9 
:  5 

!   7 

:  28 

i       i       i 

!9. 7752.41! 
! 10. 4!. 2951 
! 25. 4! 4. 99! 
!18.6!6.8  ! 

704 
695 
700 
690 

11 
36 

62.5 
92.5 

13 
18 
18 

6 

i      ■ 
i 8. 45! 1.36 
! 12.6! 1.57 
I21.6S9.15 
! 20. 7! 9. 27 

Depth  of  water  column  changed  from  93cm  to  100cm  to 
shift  second  resonance  from  593  to  544  Hz,  and  to 
shift  third  resonance  from  752  to  695  Hz. 


If  all  of  the   measurements   taken   are  grouped  together 

neglecting  frequency  and   coil   position   in  the  calibrator, 

the  blind  average  of  the  sensitivities  are: 

-3 
Coil  #1    M   =  13.8  X  10    I/rad//iPa 
CI 


with  a  standard  deviation  of: 


-3 


0~   =   7.5   X    10        /irad//iPa 
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and 


-3 
Coil  #2    M   =  13.9  X  10    U  rad//i  Pa 
C2 

with  a  standard  deviation  o-fi 

-3 
(7  =  7.7  X  10   fj,  rad//iPa 

As  can  be  seen  -from  Table  III  and  the  blind  averages,  the 
individual  coil  sensitivities  are  matched  well  within 
experimental  error. 

Table  III  indicates  an  increase  in  coil  sensitivity  as 
depth  is  increased  with  exception  of  very  deep  depths  where 
proximity  of  the  sensors  to  the  source  appears  to  have 
effected  the  measurements. 

C.   GRADIENT  SENSOR  SENSITIVITY 

A  gradient  sensor  was  formed  by  placing  the  two  fiber 
coils  in  the  calibrator  separated  by  a  distance  of  10  cm. 
The  LC-10  hydrophone  is  placed  in  the  center  of  the  coils 
approximately  5  cm  from  each  coil.  The  LC-10  hydrophone  and 
both  of  the  coils  are  then  submerged  in  water  and  exposed  to 
an  acoustic  field  and  an  established  acoustic  pressure 
gradient.  If  the  pressure  sensitivities  of  the  two  coils 
are  closely  matched,  the  pair  does  not  respond  to  the 
pressure,  but  to  the  spatial  pressure  gradient.  Such  a 
gradient  sensor  can  be  utilized  to  detect   the  direction   of 
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an  incident  acoustic   wave   as   discussed  under  conventional 

gradient  hydrophones  in  chapter  II. 

The  Jll   transducer  source  of  the  acoustic  calibrator  is 

driven  at  a  frequency  to  produce   a   well   defined   standing 

wave.   The  LC-10  hydrophone/  dual  coil  setup  is  submerged  to 

a  location  of  a  pressure   minimum,   defined  as  a  point  where 

the  voltage  from  the  LC-10  hydrophone  is  at  a  minimum.   With 

the  LC-10  at  this  position  the   two   coils   are   located   on 

either  side  of  the   pressure   minimum   and   therefore   in   a 

region   of   maximum    pressure    gradient.     The    external 

instrumentation   arrangement  is  the  same  as  when   individual 

coil   sensitivities  were  obtained.   The  driver  amplitude   is 

increased   until  the  amplitude  of  the  fundamental   frequency 

reaches  its  first   zero   point.   However,  the  drive  voltage, 

not   the   LC-10   voltage   was   recorded.    Since   the   LC-10 

hydrophone  is  at  a  minimum   pressure   position,   its   output 

voltage  is  quite  low  and,  therefore,   the   pressure   in   the 

region  of  the  coils  could  not   be   directly   computed   using 

just  LC-10  sensitivity. 

However,  by  knowing  LC-10  voltage   (V  )   at   a   pressure 

H 
maximum  for  a  set  drive  voltage  a  correlation  between   LC-10 

voltage  and  Jll  drive   voltage   could   be  calculated.   Using 

LC-10  sensitivity  (M  ) ,  a  value  of  pressure   per   Jll   drive 

H 
volts  is  obtained.   This  is  shown  in  the  following  formula: 


P  /Jllvolts  =  CV  at  P    D/CM  *JllvoltageD 
o  H     max     H  (4.3) 
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Recall  -From  Chapter  II,  equation  (2.6),  the  pressure  at 
a  distance  x  above  and  below  <+i")  a  pressure  minimum 
(assuming  x  <<A>  is: 


P   =  +P  kx 
+     o 


(4.4) 


Replacing   P    of   equation   (4.4)  with  equation   (4.3)   and 

o 
multiplying  by  the  drive  voltage  needed  to  reach   the   first 

zero  of  the  fundamental  (V    )  the  pressure  at   a  distance  x 

JU 
isa 


P   ■    C(P   /Jllvolts)*V         3kx 
±  o  JU 


(4.5) 


With  the  coils  equally  spaced  on  either  side  of  the 
pressure  minimum,  the  difference  in  pressure  is  twice  the 
pressure  at  one  of  the  coil  locations.  The  pressure 
gradient  is  the  pressure  difference  divided  by  the 
separation  between  the  coils,  and  is  given  by  the  following 
formula: 


Vp  =  2P  /^x 


(4.6) 
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Gradient  sensitivity  is  calculated  directly  by  the 
-following  -formula,  when  the  pressure  gradient  is  su-f-ficient 
to  null    the  amplitude   of    the   -fundamental: 

M   -   3.83rads/Y7p 
D 

(4.7) 

The  units  o-f  gradient  sensitivity  are  radians/  fj.   Pa/cm. 

Gradient  sensitivity  can  also  be  calculated  indirectly 
by  using  single  coil  sensitivity  values.  Indirect  gradient 
sensitivity  can  be  expressed  as  -followsi 

(4.8) 

where  <T)   and  (p     are   the  phase  shifts  in  coil  1  and  coil  2 

1        2 
respectively.   Using   the   pressure   value   as   expressed  in 

equation  (4.5)  the  phase  shifts  may  be  expressed  as: 

<b      -  P   *  M 
1     +     CI 


(4.9) 


and 


(b     =  P   *  M 
2     -     C2 


(4.  10) 
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where  P    and   P    are   the  pressure  amplitudes  at  the  sensor 

+         - 
coil  locations  and  M    and  M    are  single  coil  sensitivities 

CI       C2 
obtained   from   Table   III.  Inserting   equations   (4.9)   and 

equation  (4.10)  into  equation  (4.8)  and  using  the  fact   that 

P   =  -P   when  the  sensor  is  at  a  pressure  node,  the  indirect 

+      - 
gradient  sensitivity  can  be  expressed  as: 

M    *  P  CM    +  M   ]/Up 
ID    ±   CI     C2 

(4.11) 

Substituting  equation  (4.6)  into   equation  (4.11)  results  in 
an  indirect  gradient  sensitivity  formula  as  follows: 

M    =  CM    +  M   D/C2/^x3 
ID      CI     C2 

(4.12) 

where  l\  X  is  the  separation   between   the  two  coils  and  this 
distance  is  assumed  to  be  smaller  than  a  wavelength. 

Table  IV  is  a  summary  of  direct  and  indirect  gradient 
sensitivity  measurements  calculated  from  dual  coil  and 
single  coil  data  respectively.  Measurement  of  direct 
gradient  sensitivity  was  calculated  at  three  different 
calibrator  resonances  and  different  depths  associated  with 
pressure  minima  (maximum  pressure  gradient).  Indirect 
gradient  sensitivity  is  also  calculated  at  three  different 
calibrator   resonances   but   uses   the   average   sensitivity 
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obtained  at  a  pressure  maxima  -For  the  single  coils  divided 
by  the  desired  distance  between  the  coils  (equation  (4.12)). 
Error  analysis  for  direct  gradient  sensitivity  can  be 
made  by  taking  the  square  root  o-f  the  sum  at  squares  o-F  the 
constituents  o-F  equation  (4.6).  This  may  be  expressed  as 
follows:  // 


M 


2 


L  JIIV 


(4.13) 


where  C  is  the  experimental  average  speed  of  sound  in  the 
calibrator  as  shown  in  Table  I.  All  other  symbols  are  as 
explained  in  equation  (4.3)  and  equation  (4.5). 

Error  analysis  for  indirect   gradient  sensitivity  may  be 
expressed  as  follows: 


'dM 


M 


'ID 


I 


2   -1/2 


(4.14) 


where  <M>  is  the  average  deviation  for  single  coil 
sensitivities  and  M  is  the  average  of  the  single  coil 
sensitivities  as  presented  in  Table  III. 

D.   ANALYSIS 

A  general  analysis  of  single  coil  sensitivity  compared 
to  published  data  CRef.63  can  be  shown.  From  Table  III 
typical  sensitivities  Arm      on   the   order   of  approximately 
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-2 
10   Zirad///Pa  per   hydrophone  for  10  meters  of  fiber.   This 

results  in  a  sensitivity  of  10   /I  rad//i  Pa/m.  For  an  optical 

wavelength   of   0.83   microns,   the  experimental   value  of 

A  0/0  A P  is  approximately   130   x  10     pascal   .   For  the 

same  diameter  fiber  (0.4mm),  the  calculated  value  is  greater 

than  published  values  CRef.  53  by  a   factor  of  four  to  five. 

This  variation  could  be  explained   in   terms   of   the   large 

variations   in   observed   physical   properties    of    Hytrel 

coatings  resulting  in  a  more  sensitive  fiber. 

As  shown  in  Table  III,  the  single  coil  sensitivities  of 
coil  1  and  coil  2  at  the  same  depth  appear  to  be  closely 
matched  within  experimental  error.  The  sensitivity  of  the 
individual  coils  are  not  matched  over  a  depth  range  in  the 
acoustic  calibrator.  The  data  indicates  that  the 
sensitivity  of  the  single  coil  is  a  function  of  depth  with 
possible  near  field  effect  occuring  close  to  the  source 
causing  an  anomaly.  No  reason  could  be  found  for  this 
result  but  follow-on  experiments  will  investigate  this  depth 
dependent  sensitivity  of  single  coils. 

Measured  pressure  gradient  sensitivity  and  calculated 
pressure  gradient  sensitivity  as  shown  in  Table  IV  in  the 
previous  section  are  in  agreement  within  experimental 
error.  The  measured  pressure  gradient  sensitivity  obtained 
from  data  at  pressure  minima  is  shown  to  be  more  constant 
with  frequency  and  depth  than  the  calculated  pressure 
gradient  sensitivity  based   on   single   coil  measurements  at 
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pressure  maxima.  Where  the  single  coil  data  at  Table  III  is 
a  -function  of  depth  the  measured  sensitivity  of  the  dual 
coil  setup  appears  to  be  relatively  independent  of  depth 
with  an  average  value  of  0.124+0.015  /J.  rad//i  Pa/cm. 

Final  analysis  of  the  pressure  gradient  hydrophone  was 
conducted  by  connecting  the  output  of  the  photodiodes  to  the 
phase-lock  loop  locking  circuit.  The  dual  fiber  coils  with 
a  LC-10  hydrophone  in  the  center  was  lowered  in  the  acoustic 
calibrator.  The  voltage  from  the  locking  circuit  and  the 
LC-10  hydrophone  were  recorded.  Figures  4.6  and  4.7  are 
plots  of  the  recorded  data  at  593  Hz  and  752  Hz 
respectively.  As  shown  in  Figures  4.6  and  4.7  as  LC-10 
hydrophone  proceeds  to  a  minimum  gradient  hydrophone  voltage 
proceeds  to  a  maximum  and  as  LC-10  hydrophones  proceeds  to  a 
maximum  gradient  hydrophone  voltage  proceed  to  a  minimum. 
The  data  to  produce  Figures  4.6  and  4.7  is  in  good  agreement 
with  theory  except  for  data  recorded  near  the  source  or  near 
the  air — water  boundary. 
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Figure   4.6        Demodulated   Gradient   Hydrophone 
Output    at    593   Hz. 
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Figure  4.7    Demodulated  Gradient  Hydrophone 
Output  at  752  Hz. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

Since  their  invention  in  the  late  1970's  CRe-f.  3  and  4], 
many  applications  -for  fiber  optic  sensors  have  been 
suggested.  In  the  ArsA  of  low  frequency  (audio)  sound 
transduction  in  the  sea  (hydrophones) ,  it  was  felt  the  most 
reasonable  application  of  this  class  of  sensor  was  as  a 
dipole  sensor.  The  potentially  high  sensitivity  could  be 
exploited  since  the  output  of  a  dipole  sensor  is  the  small 
difference  of  two  relatively  large  signals.  Secondly,  an 
interf erometric  (homodyne)  detector  is  intrinsically 
differential  since,  in  this  case  the  interference  is  created 
by  the  difference  of  the  optical  path  lengths  in  the  two 
arms  of  the  interferometer.  An  additional  benefit  of  this 
differential  operation  is  the  high  degree  of  common  mode 
noise  rejection  which  maybe  gained  by  having  the  "sensor" 
and  "reference"  arm  in  the  same  environment  until  the  very 
end  of  the  sensor. 

The  laboratory  investigations  outlined  in  this  report 
have  demonstrated  the  feasibility  of  this  application. 
Because  the  fiber  optic  system  is  essentially  a  phase 
modulated  rather  than  amplitude  modulated  system  (equation 
2.22),  it  is  difficult  to  make  an  absolute  comparison 
between  it  and  a  piezoelectric  pressure  gradient  sensor. 
Before   the   impedance   converting,   20   dB   gain,   internal 
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amplifier,  the   USRD   type   G63   standard   pressure  gradient 

-9 
hydrophone  CRef.  133  has  a  sensitivity  of  7.5X10   V/^iPa/cm. 

The  -Fiber  optic   system   described   here  using  piezoelectric 

phase  shifter  (Chapter  III  section   A-4)   and   the   feedback 

locking  circuit  (Chapter  III  section  D-8)  had  an   output   of 

305+8  mvYradian  (Chapter  IV  section   A   and   Fig.  4.3).  This 

combined    with    an    average    gradient    sensitivity    of 

0.124+0.015  U rad/U Pa/cm   gives  an  overall  sensitivity   for 

""8 
the   fiber  system  of  3.8+0.5X10   V//i  Pa/cm.  By  this   method, 

the  fiber  system  has  a  14  dB  advantage   over   the   USRD  type 

G63. 

At  this  stage,  further  work  in  this  area   is   justified. 

The  directions  that  should  be  pursued  are  the  ruggedizing  of 

the  sensor,  and  the  development   of   a   package   which   will 

incorporate  the   sensor  and  couplers  in  a  single  submersible 

unit,  possibly  to  be  followed   by   the   integration   of   the 

laser   and   demodulators   in   the   same   unit.    Alternative 

demodulation  schemes  which   do   not   require  a  piezoelectric 

fiber  stretcher  to  maintain  "quadrature  lock"  should  also  be 

investigated.    Finally,    direct    comparison    with   other 

directional  dipole   sensors  such  as  those  on  DIFAR  sonobuoys 

should  be  made  under  real  sea  conditions. 
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APPENDIX  A 

DATA  ACQUISITION  PROGRAM 

This  program  was  written  to  use  a  HP-85F  computer  to 
gather  data  more  accurately  and  quickly  than  by  manual 
means.  This  program  was  used  a-fter  initial  construction  and 
prior  to  sensitivity  experiments  to  obtain  data  -for 
computation  o-f  interferometer  characteristics  (see  Figure 
4.1).  Chapter  IV  Section  A  gives  a  detailed  explanation  o-f 
how  this  program  was  utilized  for  computerized  data 
acquisition. 
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10  OPTION  BASE  1 

20  SHORT  F(3),  A(250),  B(250,3),  C(250,3>,  E(250) 

30  REM  »*  OBTAIN  INITIAL  VALUES  FOR  PARAMETERS  ** 

40  DISP  "ENTER  FUND  FREQ" 

50  INPUT  F(l) 

60  DISP  "ENTER  SPECT  ANAL  START  FREQ" 

70  INPUT  Fl 

80  DISP  "ENTER  START  AMPLITUDE  IN  MVOLTS" 

90  INPUT  A(l) 
100  DISP  "ENTER  END  AMPLITUDE  IN  MVOLTS" 
110  INPUT  VI 

120  DISP  "ENTER  INCREMENT  VALUE  IN  MVOLTS" 
130  INPUT  I 
140  F(2)«2*F(1) 
150  F(3)«3*F(1) 

160  REM  ♦♦  PRESET  AND  INITIALIZE  HP-3314  AND  HP-3582  ** 
170  OUTPUT  707  ; "PR" 

180  OUTPUT  707  ; "FR" ,F ( 1 ) , "HZAP,A( 1) , "MVFU1 " 
190  OUTPUT  711  ;"PRS" 

200  OUTPUT  711  ; " AS8SL 1 MN 1 MD4SP 1 OAD " , F 1 , " SC 1 A V4NU4 AM2 " 
210  J=(V1-A(1))/I+1 

220  REM  ♦♦  DATA  ACQUISITION  FOR  SET  DRIVE  AMPLITUDE  ** 
230  FOR  R=l  TO  J 
240  OUTPUT  711  ; "RE" 
250  WAIT  11500 
260  ENTER  723  ;  V 
270  E(R)+V 

280  REM  **  STORE  AMPLITUDE  OF  FUND  AND  1ST  TWO  HARMONICS  ** 
290  FOR  N=l  TO  3 
300  M=(F(N)-F1)/ 1000*250 
310  OUTPUT  711  ;"MP" ,M 
320  OUTPUT  711  ; "LMK" 
330  ENTER  711  ;  B  (R,N)  ,C  (R,N) 
340  NEXT  N 
350  D=R+1 
360  A(D)=A(R)+I 

370  OUTPUT  707  ; "AP" ,A(D) , "MV" 
380  NEXT  R 

390  REM  **  OUTPUT  DATA  TO  PRINTER  ** 
400  PRINTER  IS  701,80 
410  PRINT  "FUND  FREQ  (HZ)  =",F(1) 
420  PRINT  'SOURCE  VOLTS      LC10" 
430  PRINT  "MILLIVOLTS (P-P)  VOLTS (RMS)  VF1 (RMS)  VF2(RMS) 

VF3(RMS)" 
440  FOR  R=l  TO  J 

450  IMAGE  50,12X,D.DDE,4X,D.DDE,4X,D.DDE,4X,D.DDE 
460  PRINT  USING  400  ;  A (R) ,E (R) ,B (R, 1 ) ,B(R,2) ,B (R,3) 


97 


470  NEXT  R 

480  REM  *♦  RESET  INSTRUMENTATION  TO  PRESET  CONDITION  ** 

490  OUTPUT  707  ;  "PR" 

500  OUTPUT  711  j"PRS" 

510  REM  »*  STORE  DATA  ON  TAPE  »♦ 

520  ASSIGN*  1  TO  "HYD1" 

530  PRINT#  1,1  l  A<> 

540  PRINT*  1,2  %    E<> 

550  PRINT*  1,3  |  B(,) 

560  ASSIGN*  1  TO  * 

570  REM  *♦  PLOT  DATA  ON  COMPUTER  PRINTER  ** 

580  PRINTER  IS  2 

590  FOR  K=l  TO  3 

600  GCLEAR 

610  SCALE  0,10000,0,-01 

620  XAXIS  0,1000 

630  YAXIS  0,-001 

640  PENUP 

650  FOR  1*1  TO  J 

660  PLOT  A(I) ,B(I,K) 

670  NEXT  I 

680  COPY 

690  NEXT  K 

700  GCLEAR 

710  SCALE  0,10000,0,-01 

720  XAXIS  0,1000 

730  YAXIS  0,-001 

740  FOR  K=l  TO  3 

750  PENUP 

760  FOR  1=1  TO  J 

770  PLOT  A<I) ,B(I,K) 

780  NEXT  I 

790  NEXT  K 

800  COPY 

810  DISP  "END" 

820  END 
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APPENDIX  B 


SAMPLE  RAW  DATA 


PZT  #1 


Laser  Diode 

35  mi  11 lamp 

B 

Current 

Fundamental 

-frequency  (Hz) 

372 

Source 

Fundamental 

1st  Harmonic 

2nd  Harmonic 

Voltage 

Voltage 

Voltage 

Voltage 

(volts) 

(millivolts) 

(millivolts) 

(millivolts) 

(rms) 

(rms) 

(rms) 

(rms) 

0.035 

0.783 

0.030 

0.0019 

0.053 

1.180 

0.066 

0.0039 

0.071 

1.570 

0.105 

0.0048 

0.071 

1.980 

0.182 

0.0106 

0.089 

2.320 

0.239 

0.0164 

0.106 

2.700 

0.349 

0.0241 

0.124 

3.050 

0.456 

0.0347 

0.142 

3.430 

0.544 

0.0473 

0.159 

3.800 

0.642 

0.0647 

0.177 

4.070 

0.812 

0.0830 

0.195 

4.380 

0.916 

0.1150 

0.230 

4.820 

1.  130 

0. 1460 

0.247 

5.010 

1.270 

0. 1720 

0.247 

5.340 

1.510 

0.2060 

0.265 

5.590 

1.610 

0.2540 

0.283 

5.790 

1.750 

0.2880 

0.300 

6.110 

2.010 

0.3280 

0.318 

6.310 

2.150 

0.4070 

0.336 

6.460 

2.370 

0.4520 

0.353 

6.850 

2.440 

0.5200 

0.373 

6.870 

2.460 

0.5770 

0.390 

7.  130 

2.960 

0.6490 

0.409 

7.240 

2.970 

0.7160 

0.426 

7.330 

3.320 

0.8160 

0.444 

7.370 

3.420 

0.8880 

0.461 

7.550 

3.630 

0.9650 

0.478 

7.700 

3.830 

1 . 0500 

0.497 

7.750 

3.990 

1. 1300 

99 


0.513 

7.650 

4.070 

1 . 2600 

0.532 

7.750 

4.320 

1.3100 

0.549 

7.890 

4.250 

1 . 4400 

0.567 

7.460 

3.990 

1 . 5000 

0.584 

7.730 

4.690 

1 . 6700 

0.603 

7.770 

4.360 

1 . 7400 

0.620 

7.650 

5.030 

1.8100 

0.638 

7.450 

4.720 

1 . 9400 

0.655 

7.460 

5.190 

2.0200 

0.672 

7.350 

4.910 

2.2200 

0.691 

7.280 

5.080 

2.2900 

0.707 

7.070 

5.450 

2.4000 

0.726 

6.930 

5.460 

2.5400 

0.743 

6.850 

5.650 

2.6700 

0.761 

6.590 

5.730 

2.8100 

0.778 

6.500 

5.820 

2.9500 

0.797 

6.200 

5.970 

3.0400 

0.814 

6.050 

6.050 

3. 1800 

0.832 

5.760 

6.100 

3.2900 

0.849 

5.550 

6.320 

3.3600 

0.868 

5.290 

6.230 

3.4900 

0.885 

5.070 

6.370 

3.6100 

0.902 

4.770 

6.500 

3.6900 

0.920 

4.490 

6.480 

3.8500 

0.937 

4.200 

6.540 

4.0000 

0.956 

3.870 

6.540 

4.0600 

0.973 

3.590 

6.480 

4. 1900 

0.991 

3.270 

6.480 

4.2200 

1.010 

3.000 

6.560 

4.4500 

1.030 

2.710 

6.500 

4.5600 

1.040 

2.360 

5.910 

4.6200 

1.060 

2.040 

6.210 

4.7300 

1.080 

1.750 

6.170 

4.7200 

1.100 

1.480 

6.120 

4.8400 

1.110 

1.170 

6.070 

4.8200 

1.130 

0.858 

5.550 

5.0100 

1.150 

0.617 

5.630 

5. 1200 

1.170 

0.346 

5.460 

5.0600 

1.190 

0.339 

5.320 

5. 1600 

1.200 

0.611 

5.150 

5.2500 

1.220 

0.841 

4.840 

5.2900 

1.240 

1.080 

4.710 

5.2800 

1.260 

1.390 

4.570 

5.2900 

1.270 

1.610 

4.280 

5.2800 

1.290 

1.920 

4.050 

5.2900 

1.310 

2.170 

3.630 

5.3600 

1.330 

2.370 

3.560 

5.4000 

1.340 

2.600 

3.370 

5.3400 

1.360 

2.830 

3.070 

5.2500 

1.380 

2.980 

2.840 

5.2500 

1.400 

3.260 

2.660 

5. 1900 

100 


1.420 

3.340 

2.410 

5. 1400 

1.430 

3.540 

2.210 

5.0500 

1.430 

3.620 

1.930 

4.8700 

1.470 

3.900 

1.780 

4.9600 

1.490 

3.930 

1.500 

4.8100 

1.500 

4.110 

1.310 

4.7100 

1.520 

4.200 

1.080 

4.5800 

1.540 

4.230 

0.863 

4.4800 

1.560 

4.360 

0.583 

4.3900 

1.570 

4.440 

0.366 

4.2700 

1.590 

4.510 

0.293 

4. 1200 

1.610 

4.560 

0.345 

3.9900 

1.630 

4.530 

0.601 

3.8100 

1.640 

4.410 

0.766 

3.6000 

1.660 

4.560 

1.020 

3.4900 

1.680 

4.530 

1.220 

3.3400 

1.700 

4.460 

1.420 

3. 1300 

1.720 

4.350 

1.580 

2.9800 

1.730 

4.250 

1.780 

2.7600 

1.750 

4.250 

1.930 

2.5800 

1.770 

4.160 

2.090 

2.4400 

1.790 

4.010 

2.290 

2.2300 

1.800 

3.920 

2.480 

2.0300 

1.820 

3.890 

2.720 

1 . 8400 

1.840 

3.780 

2.890 

1 . 6600 

1.860 

3.670 

3.020 

1 . 4400 

1.870 

3.530 

2.880 

1 . 2600 

1.890 

3.210 

3.370 

1 . 0400 

1.910 

3.230 

3.640 

0.8710 

1.930 

3.030 

3.630 

0.6390 

1.950 

2.820 

3.730 

0.4490 

1.960 

2.570 

3.860 

0.2430 

1.980 

2.260 

3.570 

0.0772 

2.000 

2.320 

3.990 

0. 1690 

2.030 

2.090 

3.980 

0.3640 

2.050 

1.940 

4.080 

0.5500 

2.070 

1.710 

4.010 

0.7450 

2.080 

1.510 

4.120 

0.9230 

2.100 

1.270 

4.110 

1. 1200 

2.  120 

1.090 

4.050 

1 . 2800 

2.140 

0.866 

4.030 

1 . 4200 

2.150 

0.683 

3.960 

1 . 6300 

2.170 

0.474 

3.830 

1 . 7800 

2.190 

0.320 

3.880 

1 . 9200 

2.210 

0.361 

3.820 

2.0700 

2.210 

0.524 

3.760 

2. 1900 

2.230 

0.667 

3.640 

2.3500 

2.240 

0.909 

3.370 

2.4300 

2.260 

1.080 

3.350 

2.5400 

2.280 

1.260 

3.220 

2.7100 

2.300 

1.470 

3.090 

2.7800 

101 


2.310 

1.640 

2.960 

2.9000 

2.330 

1.830 

2.810 

3.0400 

2.350 

1.990 

2.690 

3.0900 

0.370 

1.970 

2.520 

3.0700 

2.380 

2.190 

2.370 

3. 1900 

2.400 

2.430 

2.160 

3.2100 

2.420 

2.570 

2.010 

3.3000 

2.440 

2.700 

1.870 

3.3200 

2.460 

2.790 

1.740 

3.3600 

2.470 

2.900 

1.540 

3.3400 

2.490 

0.950 

1.390 

3.2900 

2.510 

3.130 

1.210 

3.3700 

2.530 

3.180 

1.030 

3.3400 

2.540 

3.180 

0.766 

3.3100 

2.560 

3.270 

0.678 

3.3100 

2.580 

3.430 

0.452 

3.3000 

2.600 

3.380 

0.320 

3. 1600 

2.610 

3.430 

0.259 

3. 1000 

2.630 

3.430 

0.253 

3.0500 

0.650 

3.450 

0.292 

3.0500 

2.670 

3.430 

0.438 

2.9700 

2.680 

3.390 

0.586 

2.8700 

2.700 

3.490 

0.800 

2.7700 

2.720 

3.470 

0.939 

2.6800 

2.740 

3.390 

1.160 

2.5200 

2.760 

3.320 

1.260 

2.4300 

2.770 

3.310 

1.510 

2.3200 

2.790 

3.190 

1.620 

2.2100 

2.810 

3.130 

1.820 

2.0500 

2.830 

3.050 

1.970 

1.9100 

2.840 

2.970 

2.130 

1 . 7800 

2.860 

2.870 

2.260 

1 . 6500 

2.880 

2.740 

2.350 

1 . 4800 

2.900 

2.620 

2.540 

1 . 3400 

2.920 

2.510 

2.600 

1. 1900 

2.930 

2.390 

2.720 

1 . 0500 

2.950 

2.270 

2.780 

0.8880 

2.970 

2.  130 

2.910 

0.7160 

2.990 

1.990 

3.000 

0.5690 

3.000 

1.790 

3.020 

0.3920 

3.020 

1.690 

3.100 

0.2460 

3.040 

1.520 

3.070 

0.0733 

3.060 

1.360 

3.070 

0.0791 

3.070 

1.200 

3.250 

0.2470 

3.090 

0.997 

3.250 

0.3600 

3.  110 

0.855 

3.220 

0.5340 

3.130 

0.686 

3.220 

0.6920 

3.140 

0.558 

3.130 

0.8300 

3.  160 

0.404 

3.100 

0.9820 

3.  180 

0.303 

3.070 

1. 1300 

3.200 

0.454 

3.020 

1 . 2800 

102 


3.210 

0.585 

3.000 

1 . 4000 

3.230 

0.748 

2.720 

1 . 5200 

3.230 

0.904 

2.860 

1 . 6400 

3.270 

1.030 

2.640 

1 . 7800 

3.290 

1.210 

2.680 

1 . 8700 

3.300 

1.330 

2.520 

1 . 9500 

3.320 

1.480 

2.420 

2.0400 

3.340 

1.550 

2.330 

2. 1500 

3.350 

1.730 

2.150 

2. 1800 

3.370 

1.870 

2.070 

2.3100 

3.390 

1.990 

1.940 

2.3500 

3.410 

2.130 

1.850 

2.4500 

3.430 

2.270 

1.710 

2.5300 

3.440 

2.330 

1.560 

2.5400 

3.460 

2.490 

1.410 

2.6400 

3.480 

2.470 

1.280 

2.6000 

3.500 

2.260 

1.130 

2.6500 

3.520 

2.650 

0.984 

2.6300 
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